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RESEARCH  PROGRAM  ON  ADDITIVES  FOR  USE  IN 
AMMONIA  FUEL  FOR  INTERNAL  COMBUSTION  ENGINES 


ABSTRACT 

A  project  was  conducted,  for  Continental  Aviation  and  Engineering 
Corporation,  to  develop  and  evaluate  additives  to  improve  the  performance  of 
the  ammonia "fueled  engines  which  the  military  proposes  to  operate  in  localities 
where  supply  of  conventional  hydrocarbon  fuel  is  difficult.  Technical  work  on 
this  project  was  active  from  March  15,  1965,  to  December  17,  1965.  This  work 
mainly  comprised  preliminary  assessment  of  the  performance  of  candidate 
■immonia -additives  using  a  bench  reactor  to  measure  tendencies  of  the  additives 
to  increase  the  oxidation  rate  of  ammonia  and/or  cause  ignition,  followed  by 
tinal  evaluation  of  these  and  ether  additives  using  CFR  single  cylinder  engines. 
Of  the  many  gaseous,  liquid,  and  solid  additives  investigated,  hydrogen  and  ace- 
tviene  were  by  far  the  best  ammonia-additives  in  the  engines ,  In  the  CFR 
.park-ignition  engine,  hydrogen  was  the  best  additive  in  terms  of  good  engine 
performance  at  low  additive  concentration  with  normal  compression  ratios,  and 
acetylene  was  the  second  best  additive.  In  the  CFR  compression-ignition  engine, 
acetylene  was  the  best  additive  in  terms  of  good  engine  performance  at  normal 
compression  ratios,  whereas  hydrogen  performed  well  at  compression  ratios 
Tomewhat  higher  than  normal. 
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RESEARCH  PROGRAM  ON  ADDITIVES  FOR  USE  IN 

msmu  smmmjmm 


The  objective  of  this  project  was  to  develop  and  evaluate  additives 
to  improve  the  performance  of  the  ammonia -fueled  spark-ignition  and  compression- 
ignition  engines  which  the  military  proposes  to  operate  in  localities  where 
supply  of  conventional  hydrocarbon  fuel  is  difficult.  Additives  were  developed 
via  a  literature  survey  to  review  current  technology  and  secure  leads,  a  study 
of  the  theoretical  mechanism  by  which  additives  might  function,  and  a  selection 
of  promising-appearing  additives  to  be  evaluated.  To  evaluate  candidate 
additives,  a  small  bench  reactor  was  used  for  preliminary  screening,  and  both 
spark  and  compression-ignition  CFR  single  cylinder  engines  were  used  to  determine 
effectiveness  of  specific  additives  on  improving  engine  performance.  This  report 
comprises  a  brief  synops<c  of  the  principal  project  activities,  plus  appendices 
which  document  the  results  in  greater  detail. 


DEVELOPMENT  OF  ADDITIVES 


Literature  Search 


This  search  relied  mainly  on  the  Chemical  Abstracts  published  by 
the  American  Chemical  Society.  Except  for  numerous  references  to  the 
catalytic  oxidation  of  ammonia  to  nitrogen  oxides,  the  literature  contained 
relatively  few  references  to  the  combustion  of  ammonia.  The  references  are 
listed  in  Appendix  A,  and  could  be  grouped  into  three  categories:  gas  phase 
kinetics  (6  references) ,  flame  properties  (6  references) ,  and  general 
combustion  in  terms  of  engineering  and  rocketry  (8  references) .  The  latter 
category  included  some  French  research  on  combustion  of  ammonia  in  motor 
vehicle  engines.  A  resulting  U.  S.  Patent  (Frejacques)  claimed  acetylene 
concentrations  in  ammonia  of  37.  (w)  to  20%,  although  10%  to  12%  evidently 
was  the  preferred  concentration.  Frejacques  noted  that  the  favorable 
solubility  characteristics  <?'.  acetylene  in  liquid  ammonia  eliminated  the 
need  for  auxiliary  contain  cc  for  additives  and  the  need  for  other  devices 
such  as  dissociators.  />,. other  U.  S,  patent  (Davis)  eovered  ammonium  nitrate 
in  liquid  ammonia  as  a  iuel  for  internal  combustion  engines. 

Study  of  Additive  BSshaniam 

The  main  problem  with  acoonia  combustion  w.  s  regarded  to  be 
relative  inertness  to  oxidation  as  compared  with  hydrocarbons.  A  successful 
ammonia-additive  should  overcome  this  inertness  by: 

1)  initiating  the  combustion  chain  at  lower  temperatures; 

2)  sustainins  the  combustion  chain  process  once  initiated. 
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l'o  provide  information  on  the  ammonia  combustion  process,  and 
thus  on  the  mechanism  by  which  an  additive  night  function,  a  bench  reactor 
apparatus  was  devised  which  would  facilitate  the  measurement  of  the  ignition 
temperature  and  level  of  oxidation  of  a  mixture  of  ammonia  and  oxygen,  and 
which  bubsequently  could  be  used  in  evaluating  amnonia -additives .  Gaseous 
ammonia  and  oxygen  were  flowed  separately  to  this  reactor,  which  was  installed 
in  a  sma.l  furnace,  and  were  heated  until  ignition  occurred.  Reactor 
temperature  vs.  elapsed  time  were  recorded  on  a  strip  chart  and  ignition 
(sudden  temperature  rise)  was  read  directly.  The  consumption  of  ammonia 
was  measured  with  an  infra-red  spectrometer.  Flow  rates  were  kept  constant 
for  all  runs,  oxygen  at  100  cc/min  and  ammonia  at  56  cc/min.  The  reactor 
and  it*  operation  are  described  in  greater  detail  in  Appendix  B,  which 
also  presents  a  schematic  drawing  of  the  reactor  and  a  photograph  of  the 
installation. 

Preliminary  work  showed  that  no  ignition  could  occur  up  to  700*C 
at  atmospheric  pressure  with  any  ratio  of  ammonia  to  oxygen:  what  did  occur 
was  considerable  oxidation.  The  oxidation  rates  expressed  in  terms  of 
ammonia  consumed  were  determined  at  three  levels  of  ammonia  concentrations 
and  over  the  temperature  range  500-800*0.  The  results  are  illustrated  by 
Fig.  1.  The  rate  of  oxidation  is  seen  to  be  inversely  proportional  to  the 
ammonia  concentration.  Moreover,  at  high  concentrations,  the  oxidation 
goes  through  a  region  cf  negative  temperature  dependence. 

Another  experiment  was  made  to  determine  the  role  of  oxygen  in 
the  combustion  process .  Experiments  were  performed  in  the  presence  and 
absence  of  oxygen.  It  was  found  that  the  rate  of  ammonia  consumed  in  the 
presence  of  oxygen  was  7  times  greater  than  the  rate  in  the  absence  of 
oxygen.  This  indicated  that  oxygen  does  not  only  react  with  the  hydrogen 
produced  by  thermal  Initiation,  but  plays  an  important  role  in  initiating 
chain  reactions  of  ammonia . 

Selection  of  Additives 

In  light  of  the  preceeding  information,  additives  were  selected 
according  to  possible  use  with  ammonia  in  spark  ignition  engines  and  in 
compression  ignition  engines.  For  spark  ignition  engines,  the  primary 
interest  was  in  gaseous  additives  which  could  be  introduced  with  gaseous 
ammonia,  although  some  liquid  additives  were  considered  for  introduction 
via  Lae  air-ammonia  stream.  This  selection  thus  mainly  comprised  gases 
which  could  be  prepared  on-site  in  any  locality,  sor.1  as  hydrogen,  and  gases 
which  could  be  prepared  on-site  in  some  localities >  as  butane. 
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For  compression  ignition  engines,  the  primary  interest  was  in 
additives  that  might  be  effective  at  low  concentrations  and  that  would  be 
soluble  in  liquid  ammonia,  although  gaseous  additives  also  were  considered 
for  introduction  via  the  intake  air.  This  selection  was  mainly  bas'.-d  on 
providing  a  good  representation  from  the  classes  of  compounds  which  had 
at  some  time  been  claimed  as  cetane  improvers  as  well  as  providing  repre¬ 
sentation  between  compoun'  -  which  might  be  classed  as  "catalytic  oxidizers" 
and  as  "explosives". 


EVALUATION  OF  ADDITIVES  IN  BENCH  REACTOR 
Additives  For  Spark  Ignition  Engines 

Eight  gases  and  two  liquids  were  tested.  Two  of  the  gases, 
hydrogen  and  ozone,  appeared  to  have  a  strong  effect  on  the  oxidation  of 
ammonia,  whereas  the  remainder  had  either  a  slight  or  no  effect.  Both 
liquids,  a  mixture  of  iso-octane  and  n-heptane,  and  a  light  lubricating 
oil,  had  strong  positive  effects  on  oxidation.  Overall,  the  scope  of  the 
work  on  additives  for  spark  ignition  engines  burning  ammonia  was  quite 
Limited  because  the  favorable  performance  of  hydrogen  and  acetylene  was 
well  known ,  whereas  emphasis  on  additives  for  compression  ignition  engines 
was  desired  by  Continental.  Detailed  results  of  the  bench  reactor  work  on 
additives  for  spark  ignition  engines  are  given  in  Appendix  C. 

Additives  For  Compression  Ign.'tion  Engines 

A  total  of  45  additives  were  tested.  All  were  considered  to  be 
-.oluble  in  liquid  ammonia,  and  several  were  tested  at  several  concentrations. 
Of  this  total,  one,  styphnic  acid,  was  found  to  indu  »  ignition  at  concentra¬ 
tions  down  to  0.05%(w)  of  ammonia.  Six  additives  at  0.5%  substantially 
reduced  ignition  temperature  as  did  6  additives  at  1.07.,  2  at  1.57.,  and  4 
at  5.0%.  Although  not  inducing  ignition,  11  additives  at  5.0%  or  less 
significantly  increased  the  oxidation  of  ammonia.  Detailed  results  of 
this  work  are  given  in  Appendix  D. 


EVALUATION  OF  ADDITIVES  IN  ENGINE 

The  CFR  single -cylinder  test  engines  and  installations  are  described 
in  Appendix  E,  along  with  an  account  of  mechanical  problems  encountered  which 
were  believed  to  be  peculiar  to  operation  on  ammonia  and  additives .  The  per¬ 
formance  of  the  various  additives  in  the  engine  can  be  summarized  as  follows: 

Spark  Ignition 

Ten  gases  and  fourteen  liquids  were  tested  as  ammonia -additives  in 
the  spark  ignition  CFR  engine  at  900  and  1800  RPM  and  8.0  to  1  compression 
ratio,  using  35°  spark  advance  and  best  power  air-fuel  ratio.  The  results 
are  shown  in  Table  I.  Engine  performance  on  ammonia  only  and  on  hydrogen 
anly  was  checked  first  for  reference  purposes;  combustion  was  achieved  at 
)00  RPM  but  not  at  1800  RPM.  Of  the  gases,  hydrogen  was  by  far  the  best 
idditive,  permitting  fairly  good  engine  operation  at  900  and  1800  RPM  at 
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concentrations  of  less  than  1.5%(v).  Acetylene  was  the  next  best  additive, 
requiring  about  6%  to  give  good  engine  operation  at  both  speeds.  Ethylene, 
ethane,  and  butane  permitted  good  performance  at  relatively  low  concentrations 
at  900  RPM,  but  concentration  had  to  be  increased  substantially  to  obtain 
operation  at  1800  R7M.  The  remaining  gases,  nitrous  oxide,  nitrogen  dioxide, 
carbon  monoxide,  Freon-12,  and  ozone,  were  ineffective.  The  liquid  additives, 
although  several  permitted  good  engine  operation  at  both  900  and  1800  RPM,  all 
required  use  at  unreasonably  high  (over  107.)  concentrations. 

Hydrogen  also  was  checked  as  an  additive  at  several  concentrations 
at  varying  spark  advance  and  compression  ratio.  The  results  for  900  RFM  are 
shown  in  Table  II,  and  for  1800  RPM  are  shown  in  Table  III.  Hydrogen  was  an 
effective  ammonia -additive  in  this  engine  at  concentrations  of  less  than  1.0X, 
provided  compression  ratio  was  above  12:1  and  spark  advance  was  more  than  50®. 
Checks  on  carbon  monoxide  and  nitrous  oxide  were  made  also;  fairly  good  engine 
operation  could  be  obtained,  even  at  1800  RPM,  provided  additive  concentration, 
compression  ratio  and  spark  advance  were  sufficiently  high. 

Compression  Ignition 

This  work  comprised  an  evaluation  of  the  performance  of  gases  added 
to  the  intake  air,  and  of  several  liquid-ammonia  additives  selected  from  the 
bench  reactor  program.  The  work  on  the  inducted  gases  was  conducted  first,  to 
simultaneously  evaluate  the  gases  as  ammonia  additives  and  pinpoint  the  best 
engine  operating  conditions  to  use  for  evaluating  the  liquid -ammonia  additives. 

Of  gases  added  to  the  intake  air  while  ammonia  was  being  injected, 
hydrogen,  normal  butane,  and  acetylene  all  facilitated  combustion  at  some 
engine  operating  condition.  The  gases  generally  were  tested  at  20:1,  25:1 
and  30:1  compression  ratio,  and  at  900  and  1800  RPM.  The  results  are  suuna- 
rized  in  Table  IV.  Combustion  with  ammonia  and  hydrogen  was  achieved  at  25:1 
and  at  30:1  compression  ratio  at  both  900  and  1800  RPM;  performance  improved 
as  either  compression  '.atio  or  speed  was  increased.  Combustion  with  ammonia 
and  r.?nnal  butane  was  achieved  at  all  compression  ratios  and  speeds  (not 
tested  at  30:1  CR  and  1800  RPM)  but  butane  consumption  was  excessive  at*. 

1800  RPM.  Combustion  with  ammonia  and  acetylene  was  achieved  at  20:1  and 
25:1  CR  and  900  and  1800  RPM  (not  tested  at  30:1  CR) .  Performance  improved 
as  either  compression  ratio  or  speed  was  increased. 

The  additives  selected  from  the  bench  reactor  program  were  engine 
tested  at  25:1  CR  and  900  RPM,  the  condition  where  all  inducted  gases  per¬ 
formed  most  similarly.  Of  the  "non-explosives",  ammonium  nitrate  was  tested 
at  concentrations  of  1.0  and  5.07»  of  r.n  additive-liquid  anmonia  mixture,  and 
ammonium  perchlorate  was  tested  at  1.0%.  Styphnic  acid  was  selected  as  the 
most  effective  "explosive"  and  v-.a  tested  at  0.1%.  During  the  testing  of  these 
additives,  the  engine  was  closely  observed  for  any  evidence  of  the  additive 
affecting  combustion,  and  was  checked  for  fireability  by  bracketing  runs  on 
the  additives  with  runs  using  inducted  acetylene.  All  of  these  liquid -ammonia 
additive3  proved  to  be  completely  without  discernible  tendency  to  initiate 
combustion  of  suuuon la  iu  engine. 
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CONCLUSIONS 

1)  For  a  CFR  single-cylinder  spark  ignition  engine  burning 
ammonia ,  hydrogen  was  the  best  additive  in  terms  of  good  engine  performance 
at  low  additive  concentration  with  normal  compression  ratios,  and  acetylene 
was  the  second  best  additive. 

2)  For  a  CFR  single-cylinder  compression  ignition  engine  burning 
ammonia,  acetylene  was  the  best  additive  in  terns  of  good  engine  perfozuence 
at  normal  compression  ratios,  whereas  hydrogen  performed  well  at  compresaion 
ratios  somewhat  higher  than  normal. 
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PERFORMANCE  OF  AMiONIA  AND  ADDITIVES  IN  SPARK  IGNITION  ENGINE  AT  «00  ini 
_ AT  VARYING  SPARK  ADVANCE  ANT*  SSf  *T  9°°  ** 

r  ii«»i  ■muirtiii  —a 

Horaepower  at  Beat  Air/Fuel  Ratio<l)  for  the 
Following  Compression  Ratios  and  Spark  Advance* 


4.5:1 


COMPRESSION  RATIO 


12.3:1 


Ammonia 


Ammonia  +0.6%  Hydrogen 


Ammonia  +  1.1%  Hydrogen 


Ammonia  +  2.3%  Hydrogen 


Ammonia  +  3 ,4%  Hydrogen 


SPARK 

advance 


Hydrogen 

4  2.4  2.0 

7:1  8.5:1 

Ammonia  +  59%  Nitroua  Oxide 
5.1 

Ammonia  +41%  Carbon  Monoxide 


(1)  7.7 

5.3 
7.9 

CJD:bjb 

1/4/66 


:i  for  ammonia-hydrogen 

•  1  for  ftWNMn^  .• — A  Ja.  «. — s 

‘  - - - UAAUC j  aUU 

:1  for  ammonia-carbon  monoxide. 
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PERFORMANCE  OP  AWOKIA  AND  ADDITIVES  IN  STARK  IGNITION  ENGINE  AT  1800  MM 
_ AT  VARYING  SPARK  ADVANCE  AND  COMPRESSION  RATIO  . . - 

Horsepower  At  Best  Air /Fuel  Ratios 
For  The  Following  Compression  Ratios 
And  Saark  Advances  .  _ 


COMPRESSION  RATIO 


6 *  '• 


Ammonia  + 

.7%  H2 

0 

0 

0 

0 

7.2  (55) 

7.3  (50) 

0 

1.3%  H2 

0 

0 

6.7  (50) 

6.9  (40) 

7.4  (45) 

7.6  (40) 

11 

2.6%  H2 

0 

6.4  (60) 

7.0  (35) 

7.5  (35) 

7.7  (30) 

8.1  (30) 

3.8%  H2 

3.4  (45) 

6.7  (40) 

7.5  (35) 

7.8  (30) 

8.1  (30) 

Very 

i] 

2.6%  H2 

m 

6.4  (40) 

7.0  (35) 

7.5  (35) 

0 

erratic 

0 

a 

19%  CO 

0 

0 

0 

6.3  (35) 

8.4  (35) 

8.3  (35) 

15%  N^O 

0 

0 

0 

0 

4.7  (45) 

7.5  (45) 

.0 

49%  NjO 

0 

0 

0 

0 

0 

6.7  (35) 

i) 

.3%  03 

0 

0 

0 

0 

0 

0 

1) 

0 

Note:  Numbers  in  parenthesis  are  spark  advances  in  cBTDC. 

Unless  otherwise  indicated,  the  spark  plug  gap  is  .035  inches.  n 


v^)  7.9:1  for  ammonia -hydrogen,  7.2:1  for  ammonia-carbon  monoxide,  and 

7.6:1  for  ammonia-nitrous  oxide. 


.  i/66 


i 

I 


1/4/66 
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DESCRIPTION  OF  BENCH  REACTOR  AND  OPERATION 


The  apparatus  shown  in  Figs.  B-l  and  B-2  consists  of  a  reactor, 
a  furnace  for  heating  the  reactor  and  preheating  the  oxygen,  and  an  additive 
solution  atomizer.  Ammonia  and  oxygen  entered  the  apparatus  separately.  The 
ammonia  entered  at  the  top  of  the  atomizer  at  a  rate  of  56  cc/min.  After 
mixing,  the  oxygen  and  ammonia  and  additive,  (sdien  used)  passed  into  the 
furnace  tfiich  was  programmed  for  a  linear  temperature  increase.  The  ammonia - 
oxygen  mixture  flowing  into  the  reactor  was  heated  progressively  to  higher 
temperatures  until  ignition  occurred.  The  temperature  of  the  oxidation 
reaction  was  monitored  by  a  thermocouple  placed  in  the  reactor.  Temperatures 
were  measured  in  the  reactor  (T^)  indicated  by  a  "kick"  on  the  temperature 
recorder.  Ihe  effluent  gas  was  passed  directly  into  a  Perkin  Elmer  Model  112 
Infra  Red  Spectrometer  which  measured  the  concentration  of  ammonia. 

When  testing  additives,  the  additive  was  fed  to  the  atomizer  by  a 
syringe  pump,  generally  at  a  rate  of  0.01  cc/min.  The  additive  solution 
entered  the  apparatus  through  a  1/16"  O.D.  stainless  steel  tube  with  a  sharp 
pointed  tip  which  was  maintained  at  a  potential  of  about  15  KV.  The  solution 
was  charged  as  it  came  to  the  tip  and  atomized  as  it  was  pulled  into  the 
high  potential  field.  Atomization  was  achieved  by  charging  the  additive 
solution  to  a  high  potential  and  exposing  it  to  a  high  potential  field. 

Additive  solutions  were  made  up  previously,  at  the  required 
concentrations.  In  most  cases  the  solvent  used  was  water.  Ammonium  hydroxide 
solution  was  used  with  hard  to  dissolve  compounds.  In  the  case  of  insoluble 
liquids,  a  detergent  was  added  and  an  emulsion  was  used. 

With  additives  which  did  not  induce  ignition  below  600°C.  the  pro¬ 
grammer  was  stopped  at  600°C  and  the  system  allowed  to  equilibrate.  At  this 
temperature,  the  concentration  of  ammonia  was  observed.  This  was  compared 
with  a  similar  run  with  no  additive  present,  and  the  difference  was  ascribed 
to  the  induced  oxidation  by  the  additive. 

The  initial  work  with  gaseous  additives  was  performed  by  introduction 
of  the  gas  at  top  of  the  reactor  without  the  voltage  turned  on  the  atomizer. 
Measurement  of  oxidation  increase  was  made  by  a  differential  temperature 
method.  The  exothermicity  of  reaction  could  be  measured  by  the  difference  in 
temperature  readings  between  (Tg)  and  (T*>) . 
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APPENDIX  C 

BENCH  REACTOR  RESULTS  ON  ADDITIVES  FOR 
SPARX  IGNITION  ENGINES  BURNING  AMMONIA 


Additive 

Concentration, 

V  (v) 

Effect  on  Oxidation 
of  Amenta 

Hydrogen 

3-  5 

Strong 

Ozone 

1 

If 

n -Butane 

3 

Weak 

Butene-1 

3 

fi 

Carbon  Monoxide 

3 

19 

Chlorine 

3 

It 

Nitrous  Oxide 

1-10 

None 

Nitrogen  Dioxide 

1-10 

fi 

Iso-Octane-n-Heptane, 

80/207. 

3 

•Strong 

Lube  Oil  Base,  No.  5, 

Solvent  Extracted 

5 

u 
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APPENDIX  D 
TABLE  D-I 

ADSmjgS„FOR,LI<yiID  AMMONIA  WHICH  INDUCED-  IGNITION  IN  THE  BENCH  BEAETOft 


No  Additive 


additives  at  <0.5%  (w)  of  aranonia: 
Styphnic  Acid  (0.35%) 

"  ”  (0.25%) 

”  "  (0.10%) 

"  "  (0.05%) 

Additives  at  0.5%  (w)  of  ammonia: 
Styphnic  Acid 
Trinitroraetacreaol 
Nitroguanidine 
Trinitrotoluene 
Picric  Acid 
Ammonium  Picrate 

Additives  at  1.0%  (w)  of  ammonia: 
Ammonium  Perchlorate 
Sodium  Azide 
Trinitrobenzene 
Trinitromethane 
Isoamyl  Nitrite 
Amyl  Nitrate 

Additives  at  1.5%  (w)  of  anoonia: 
Ammonium  Nitrate 
Potassium  Azide 

Additives  at  5.0%  (w)  of  awnonia: 
Potassium  Nitrate 
Ammonium  Perchlorate 
Ammonium  Carbonate 
Sodium  Nitrite 

Additives  at  >5.0%  (w)  of  ammonia: 
Ammonium  Nitrate  at  7.5% 

Ammonium  Nitrate  at  15.0% 
Hydrczine  at  35.0% 


Ignition  Temperature, 
■■Degrees  Centigrade 

>  700 


400 

395 

520 

580 


395 

410 

430 

455 

460 

580 


480 

500 

515 

525 

570 

645 


540 

685 


450 

460 

610 

625 


535 

420 

510 
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TABLE  D-II 

ADDITIVES  FOR  LIQUID  AMHONIA 
WHICH  INCREASED  OXIDATION  IN  THE  BENCH  REACTOR 


Additives  at  1.0%  (w)  of  somoni* 
Potassium  .izide 
Aomoniura  Nitrate 

Additives  at  ?..5%  (w)  of  ammonia 
Potassium  Permanganate 
Sodium  Dichromate 
Potassium  Dichromate 
Triethylamine 

Additives  at  b  0%  (w)  of  ammonia 
Ammonium  Persulfate 
Sucrose 

Hydrogen  Peroxide 

Hydrazine 

Urea 


1 


I 

I 

I 

I 

1 

1 

! 

I 

I 

I 

I 

1 

l 

1 

l 

l 

r 

r 
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TABLE  D-III 

ADDITIVES  FOR  LIQUID  AMMONIA 
WHICH  ARE  INEFFECTIVE  IN  THE  BENCH  REACTOR 


Additives  at  0.257.  (v)  of  ammonia: 
Trinitroraetacresol 
Nltroguanidine 

Additives  at  0.51  (v)  of  ammonia: 
Trinitrobenzene 
Tetranitromethane 

Additives  at  1.07.  (v)  of  ammonia: 
Nitroethane 
Nitropropane 
Potassium  Perchlorate 
Sodium  Nitrite 
2-Nitro  2  Methyl  P;:opanoI 
2-Nitro  2  Methyl  Propanol  Nitrate 

Additives  at  2.0%  (w)  of  ammonia: 
Ethyldiazoaccca  te 

Additives  at  2.57.  (w)  of  am..  : 
Potassium  chlorate 
Ammonium  Vanadate 
Tert  Butyl  Acetate 
Tert  Butyl  Nitroso  Mercaptan 
Butyl  Mercaptan 
n-Butyl  Ether 
Acetomide 
Acetoxlme 

Additives  at  5.0%  (w)  of  ammonia: 
Ammonium  Oxalate 
Ammonium  Acetate 
Ammonium  Chloride 
Sodium  Nitrate 
Sodium  Chlorate 


i 


! 

* 


I 


1 

1 
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APPENDIX  E 

TEST  ENGINES,  INSTALLATION,  AND  MECHANICAL  PROBLEMS 
RELATED  TO  OPERATION  ON  AMMONIA  AND  ADDITIVES 

Spark  Ignition  Engine 

A  CFR  Research  method  engine  was  used.  The  standard  CFR  fuel 
system  was  omitted  and  the  engine  was  fitted  with  a  manifold  that  per¬ 
mitted  simultaneous  independent  introduction  of  ammonia  and  additive. 

Air  flow  measurements  were  made  using  a  Meriam  Laminar  Flow  Meter. 

Ammonia  flow  measurements  were  made  using  a  Fischer  and  Porter  Variable 
Area  Flowmeter.  Average  coolant  and  oil  temperatures  were  about  212*F 
and  135*F  respectively.  Mechanical  experiences  with  this  engine  were: 

1)  Excessive  leakage  of  ammonia  was  experienced  through  the  crank¬ 
shaft  and  injection  pump  drive  seals  and  the  distributor  shaft.  This 
condition  was  corrected  by  fitting  lip  seals  to  the  crankshaft,  an  "0- 
ring"  seal  to  the  distributor,  and  sealed  ball  bearing  sets  to  replace 
standard  ball  bearing  sets. 

2)  Nitrogen  dioxide  addition  resulted  in  extensive  spark  plug 
fouling  and  substantial  combustion  chamber  deposits.  The  deposits  were 
identified  as  ammonium  nitrate.  These  deposits  probably  contributed  to 
engine  wear  which  required  replacement  of  the  cylinder  and  piston.  At 
the  time  of  replacement,  the  engine  had  operated  on  ammonia  and  additives 
for  a  total  of  110  hours. 

3)  The  ignition  system  consisted  of  an  automotive  ignition  coil,  a 
4  mfd.  capacitor,  UD-16  spark  plug  gapped  at  .035",  and  a  primary  voltage 
of  12V.  No  ignition  difficulty  was  experienced  in  this  system. 

To  provide  some  bench-mark  information  on  ignition  requirements, 
oscillograms  oi  voltage  traces  from  the  high  tension  lead  wire  were  taken 
while  the  engine  was  running  with  ammonia  and  with  ammonia  plus  hydrogen. 

The  oscillograms  showed  numerous  high  voltage  spikes  following  firing,  plus 
much  hash  on  the  decay  curve,  upper  traces,  Figs.  E-l  and  E-2.  These  traces 
indicate  that  ionization  took  place  at  about  6-7  kilovolts  (1  cm  =  2  KV) . 

Compression  Ignition  Engine 

A  CFR  Cetane  method  engine  was  used  for  this  phase.  Figs.  E-3 
and  E-4.  Average  coolant  and  oiJ  temperatures  were  205-212#F  and  130-190*F 
respectively.  Air  flow  measurements  were  made  using  a  Meriam  Laminar  Flow 
Meter.  The  injection  pump  rack  was  calibrated  on  heater  oil.  Ammonia  flows 
were  then  determined  in  terms  of  equivalent  flow  on  heater  oil.  Mechanical 
experiences  with  this  engine  were: 
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1)  Piston  seizure  was  experienced  three  tines  using  cast-iron  pistons 
during  the  3hake-down  runs.  Th6  piston  seizures  evidently  were  due  to 
lack  of  adequate  lubrication  and  cooling,  and  to  difficulty  in  establishing 
proper  piston-to-cylinder  clearances  at  higher  engine  speeds,  particularly 
at  270C  *FM.  One  piston  developed  a  hairline  crack  from  the  skirt  to  the 
wrist-pin  boss.  No  further  difficulty  occurrev  when  an  aluminum  piston  was 
substituted  for  the  cast-iron  piston.  The  advantages  of  the  aluminum  piston 
apparently  lay  in:  1)  greater  clearances  used  in  fitting,  2)  greater 
strength  of  the  skirt,  3)  better  heat  conductivity  and  4)  less  risk  of 
damage  to  the  cylinder  - 

2)  Head  gasket  failure  was  a  continuing  problem.  Three  gasicets 
failed  in  rapid  succession,  Fig.  E-5.  The  failures  were  thought  to  result 
from  the  high  pressures  developed  at  elevated  compression  ratios.  Several 
attempts  to  correct  this  difficulty  were  made  by  substituting  a  homemade 
copper  gasket  for  the  composition  gasket  which  is  standard  with  the  C?R  engine. 
The  homemade  gaskets  were  not  completely  successful . 

Then  the  engine  was  modified  in  an  effort  to  prevent  further 
head  gasket  failures.  The  cylinder  was  counterbored  a  short  distance  down 
below  the  valve  deck,  and  an  aluminum  fire-ring  was  fitted  to  the  groove. 

The  ring  projected  above  the  deck  a  distance  about  equal  to  the  crush - 
thickness  of  the  standard  head  gasket.  The  fire  ring  was  cut  from  a  standard 
gasket  to  enlarge  the  hole  to  accommodate  the  aluminum  ring.  Th>*s,  the 
aluminum  ring  sealed  the  combustion  chamber  and  the  gasket  sealed  the  water 
jacket.  The  aluminum  rings  were  cut  from  sand-cast  aluminum. 

This  approach  was  reasonably  successful.  However,  the  sand-cast 
aluminum  fire  ring  eventually  failed  and  had  to  be  replaced.  The  ring 
showed  heavy  ablation  (Fig.  E-6)  in  the  same  relative  location  where  the 
head  gaskets  had  failed  previously.  A  second  ring  of  the  same  material 
was  installed.  Upon  failure,  it  was  found  that  half  of  the  ring  had  dis¬ 
appeared,  Fig.  E-7.  Several  rings  were  then  fabricated  from  2024-T351 
aluminum,  it  being  thought  that  the  higher  Young's  modulos -might  reduce 
ablation  occurring -during  periods  whin  the  cylinder  head  might  flex  away 
from  the  seal.  However,  upon  final  engine  disassembly,  a  2024-T351  ring 
was  removed  which  had  a  1-1/2  inch  segmenc  missing. 

This  suggests  that  the  gasket  failures  are  probably  not  completely 
attributable  to  pressure,  or  to  irregularities  in  the  mating  surfaces.  A 
possible  explanation  might  be  that  extremely  high  flame  temperatures  exist, 
compared  to  normal  combustion,  during  the  period  of  spontaneous  ignition  of 
the  gaseous  additives  tested.  This  theory  is  supported  by  the  evidences 
of  pilot  work  conducted  earlier  (preceeding  work  under  this  contract)  wherein 
the  engine  was  operated  at  elevated  compression  ratios  (37:1)  and  high 
jacket  (350*F)  and  inlet  air  (200*F)  temperatures  for  sustained  periods. 

No  gaseous  additives  were  used  during  this  earlier  work  and  no  head  gasket 
damage  was  experienced. 
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3)  Injector  nozzle  sticking  occurred  repeatedly  during  periods  of 
acetylene  addition.  Microscopic  examination  revealed  varnish-like  deposits 
on  the  injector  pintle.  This  apparently  was  due  to  polymerization  of 
acetylene  during  periods  when  pressures  in  the  combustion  chamber  exceeded 
the  1500  psi  Injector  nozzle  pressure  setting.  Increasing  the  injector 
pressure  to  2500  psi  reduced  the  incidence  of  nozzle  sticking.  Nozzle 
sticking  also  resulted  during  operation  on  ammonium  nitrate  and  ammonium 
perchlorate. 

4)  Two  cylinder  heads  were  cracked  during  these  tests.  Each  crack 
occurred  in  the  web  between  the  valves.  Detonation  was  often  quite  severe, 
particularly  with  hydrogen,  and  undoubtedly  caused  this  damage. 

5)  Erosion  of  the  piston  crown  area  was  observed,  plus  extensive 
nibbling  of  the  piston  corner  in  an  area  adjacent  to  the  point  where  the 
charge  enters  the  cylinder.  This  could  be  due  only  to  chemical  attack  by 
ammonia  on  the  aluminum  and/or  to  preignition. 

Examples  of  dp/dt  traces  for  heater  oil  and  ammonia  +  various 
additives  are  shown  in  Figs.  E-8  to  E-15. 

The  upper  trace  is  injector  travel  and  the  left  boundary  of  the 
oscillogram  mark*  the  beginning  of  injection.  Hence,  at  900  RPM,  we  see 
an  ignition  delay  of  17*  for  heater  oil,  37°  for  ananonia  +  n-Butane,  and 
37°  for  acetylene.  Oscillograms  for  ammonia  +  hydrogen,  at  900  RPM,  are 
not  available  due  to  failure  of  the  pressure  pick-up,  however,  the  traces 
characteristically  looked  like  those  for  ammonia  +  n-Butane  but  with  a 
high  spike  on  the  peak. 

Combustion  chamber  resonance  is  evident  with  heater  oil  and  ammonia 
4  hydrogen  at  1800  RPM.  This  contrasts  with  the  traces  for  ammonia  +  n-Butane 
or  acetylene  wherein  there  is  no  resonance. 
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FIGURE  E - 2 

AMMONIA  +  HYDROGEN 

900  RPM 
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FIGURE  E •  5 
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FIGURE  E - 8 

HEATER  OIL 

900  RPM 


FIGURE  E-9 


HEATER  OIL 

1  £  **  0  RPM 
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FIGURE  E - 1  0 

AMMONIA  +  HYDROGEN 

1800  RPM 


FIGURE  E - 1 1 

AMMONIA  +  HYDROGEN 

2700  RPM 


-  33  - 


CAE  Report  No.  1054 
Appendix  I 
Volume  II 


FIGURE  E-l  2 
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900  RPM 


FIGURE  E-1  3 
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FIGURE  E - 1  4 

AMMONIA  +  ACETYLENE 

900  RPM 


FIGURE  E.15 

AMMONIA  +  ACETYLENE 

1800  RPM 
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INTRODUCTION 

This  Research  Note  was  prepared  by  the  Allison  Division  of  General 
Motors  to  present  the  technical  results  of  ammonia  additive 
screening  tests  conducted  for  Continental  Aviation  and  Engineering 
Corporation  under  Purchase  Order  No.  RD  101033*  Continental 
Aviation  and  Engineering  Corporation  is  investigating  ammonia  as 
a  fuel  for  spark  and  compression  ignition  engines  under  Army 
Contract  No.  DA-20-113-AMC‘*05553(T) .  This  activity  is  under  the 
technical  supervision  of  the  U.S.  Army  Tank  Automotive  Center. 

The  additives  to  be  screened  in  the  Allison  flame  tubes  were  to  be 
selected  by  Continental  Aviation  and  Engineering  Corporation 
based  on  the  results  of  an  independent  investigation  by  tha 
American  Oil  Company,  Whiting,  Indiana.  Prior  to  the  receipt  of 
definitive  additive  results  from  the  American  Oil  Company’s  in¬ 
vestigations,  Allison  recommended  and  was  authorized  to  conduct 
screening  tests  with  three  oxides  of  nitrogen  and  carbon  monoxide. 
In  addition.  Continental  requested  that  screening  tests  be  con¬ 
ducted  on  ammonia  nitrate  and  hydrogen  peroxide.  The  nitrogen 
oxides  were  of  interest  due  to  their  apparent  activity  as  a  rapid 
scavenger  of  the  amidogen  radical  in  ammonia  oxidation  and 

evidence  that  NO  reacts  directly  with  NH  and  NHg  radicals  without 


prior  decomposition 


Carhon  monoxide  and  ammonium  nitrate  were 


included  based  upon  observations  from  early  studies  by  American 
Oil  (3).  Hydrogen  peroxide  was  of  interest  due  to  its  exothermic 
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decomposition  characteristics  which  would  liberate  both  energy 
and  free  oxygen  during  the  initial  combustion  process.  All  of  the 
materials  except  carbon  monoxide  were  also  potentially  producible 
in  the  field  by  the  Energy  Depot  concept. 

CONCLUSIONS: 

1.  In  5#  volume,  none  of  the  additives  tested  increased  the 
apparent  ammonia-air  flame  propagation  velocity  to  the  degree 
obtainable  by  dissociating  28#  cf  the  incoming  ammonia  fue? . 

2.  The  beneficial  effect  of  carbon  monoxide  and  the  nitrogen 
oxides  tested  are  very  dependent  upon  the  local  fuel-air  ratio 
of  the  ammonia-air  mixture. 

3.  Nitric  oxide  appears  to  offer  the  most  promise  as  a  flame 
propagation  rate  improver  over  the  fuel  equivalence  ratio 
range  of  0.9  to  1.2.  However  it  did  not  significantly  exceed 
the  effects  of  5#  acetylene  additions  at  any  mixture  ratio. 

Also  the  formation  of  solid  deposits  resulting  when  the 
additive  vapor  in  air  mixed  with  the  ammonia  fuel  could  render 
its  engine  application  difficult. 

4.  Carbon  monoxide  provided  significant  improvement  to  the 
ammonia-air  flame  propagation  rate  at  fuel  lean  (0.9  stoichio¬ 
metric)  conditions,  but  was  ineffectual  at  fuel-rich  conditions. 


-2- 


CAE  Report  1054 


TECHNICAL  DATA  REPORT  Volume  II 

~  '  ■  EDO  NO  RfcC  ’  PACE  PAGESl  REPORT  SO 


OIV. 

ALLISON 

CMC 

EDO  NO 

PAGE  PAGESl  REPOS*  SO 

.  3  oF  /  ; 

TITLE 

PREPARED 

DATE 

CHECKED 

'  APPROVED 

- 1 

5-  Flame  stability  limits  were  extended  markedly  at  fuel -rich 
(>  1.2  of  stoichiometric)  conditions  by  nitrogen  dioxide 
addition,  but  the  benefits  were  minor  in  the  stoichiometric 
region.  Undesirable  solid  deposits  were  also  experienced 
with  this  material  in  the  fuel-air  mixing  region. 

6.  No  conclusive  results  could  be  obtained  as  to  the  effect  of 
ammonium  nitrate  additions  due  to  inject?  >n  and  mixing  pro¬ 
blems  inherent  to  the  teat  burner  system. 

7.  The  burner  used  in  this  investigation  is  useful  for  screening 
the  effects  of  gaseous  or  readily  vaporizable  additives  upon 
the  propagation  rate  of  amraonia-air  flames.  However,  it  ir 
not  suitable  for  evaluating  non-volatile  additives  without 
development  effort  beyond  the  scope  of  the  present  project. 

RECOMMENDATION 

Further  additive  screening  with  the  presently  developed  test 
burner  is  recommended  only  for  readily  vaporizable  materials. 

DISCUSSION: 

Test  Eaulpment 

(4) 

The  Powling  type  flat  flame  burner  developed  for  the  previous 
ammonia  fueled  gas  turbine  burner  investigation  was  used.  The 
basic  configuration  is  shown  in  Figure  1.  The  fuel-air  mixture 
is  admitted  through  the  inlet  port  at  the  base  of  the  burner.  The 
normally  gaseous  or  vaporized  additives  are  admitted  to  the  air 
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stream  ahead  of  the  fuel-air  mixing  point.  Mixing  of  the  fuel 
and  air-additive  mixture  is  completed  as  the  mixture  progresses 
through  the  glass  bead  filled  section  in  the  base.  After  passing 
through  a  sintered  stainless  steel  flame  arrestor,  the  gases  flow 
through  the  honeycomb- glass  bead-honeycomb  section  to  flatten  the 
velocity  profile  of  the  effluent  gases.  The  burner  is  enclosed 
within  a  2.9  in.  diameter  jacket  covered  at  the  top  by  a  screen  to 
minimize  effects  of  random  air  currents  and  convective  diffusion 
of  outside  iir  into  the  reaction  zone.  The  burner  mixing  section 
and  fuel  and  air  inlet  lines  were  wrapped  with  resistance  heating 
tape  to  provide  controlled  inlet  temperatures  from  ambient  to 
500 *P.  The  air,  fuel,  and  additive  flows  were  individually  con¬ 
trolled  and  measured,  and  the  temperature  and  pressure  of  the 
incoming  mixture  was  measured  in  the  burner  mixing  section.  A  disc 
type  flat  non-diffusion  flame  adjacent  to  the  top  of  the  1/16" 
honeycomb  surface  results  at  stable  burning  conditions.  When  stable 
burning  limits  are  approached,  the  flame  normally  lifts  abruptly 
from  the  surface  and  blow-out  results.  This  condition  is  readily 
reproducible. 


The  previously  described  system  was  suitable  for  testing  the 
normally  gaseous  or  readily  vaporizable  additives.  However,  to 
evaluate  ammonium  nitrate,  modifications  were  necessary.  Electro¬ 
static  atomization  and  distribution  of  a  by  weight  01  aqueous 
ammonium  nitrate  in  the  mixing  zone  below  the  burner  flame-holding 
grid  was  attempted.  This  approach  was  selected  based  upon  recent 


s 
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work  at  California  Research  Corporation  .  The  ammonium 

nitrate  additive  supply  and  electrostatic  atomization  system  is 
schematically  shown  in  Figure  2,  and  the  final  burner  configuration 
is  presented  in  Figure  3'  A  number  of  variations  in  location  of 
the  charged  grids  and  the  additive  injector  were  explored  with  only 
limited  success  in  obtaining  the  desired  additive  solution  dis¬ 
tribution  across  th*  -eaction  zone. 

:  Due  to  the  thermal  instability  of  vaporized  hydrogen  peroxide,  the 
electrostatic  atomization  approach  was  planned  for  use  in  examining 
this  additive.  In  addition,  passivation  of  the  peroxide  inlet 
system  and  minor  modifications  to  avoid  trapping  hydrogen  peroxide 
at  any  point  in  the  system  would  also  be  necessary.  Due  to  the 
difficulties  encountered  in  effective  electrostatic  atomization  of 
the  aqueous  ammonium  nitrate,  and  the  expense  required  for  rig 
modifications  to  test  hydrogen  peroxide,  no  attempt  was  made  to 
obtain  data  on  this  material  when  information  was  received  that  the 
project  was  being  terminated  due  to  redirection  of  the  total  program. 

Experimental  Technique 

For  this  screening  effort,  all  giseous  additives  were  checked  at  a 
concentration  of  5#  by  volume-  of  fuel  flow.  Stable  burning  was 
established  and  the  fuel  a.  additive  flews  adjusted  to  the  desix*ed 
test  condition.  The  air-flov  rate  was  then  increased  or  decreased 
until  blow-out  resulted.  This  point  is  normally  a  sharp  and  j 

reprodui  ibie  condition.  An  attempt  was  mads  to  run  only  the  number 
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of  test  points  necessary  for  comparison  of  the  additive  effects.  ! 

i 

l 

Determinations  of  the  effect  of  CO,  NgO,  and  NO  additives,  and  a  , 
baseline  of  ammonia-air  were  run  at  ambient  (77*F)  inlet  temperatures 
The  NOg  (boiling  point  70*F)  additive  supply  and  inlet  line  was 
pre-heated  to  125#F  to  adequately  vaporize  the  material.  However, 
at  the  77*F  ammonia  fuel  and  air  inlet  temperature  condition, 
severe  deposition  of  solid  nitrate  deposits  occurred  at  the  fuel-air 
mixing  zone. 

Testing  at  an  inlet  air  temperature  of  300 *F  was  found  to  be  nec¬ 
essary  to  minimize  this  deposit  formation  and  to  obtain  reproducible 
results.  An  ammonia-air  baseline  was  also  obtained  at  the  300 *F  air 
inlet  temperature.  The  limited  scope  of  the  program  and  the 
indicated  additive  benefits  did  not  justify  the  additional  deter¬ 
minations  of  CO,  NgO,  and  NO  effects  at  300 *F  air  inlet  temperature 
to  provide  the  same  baseline  for  all  evaluations. 


In  the  attempts  to  evaluate  the  effects  of  the  ammonium  nitrate 
solution  additive,  the  ammonia-air  flame  was  stabilized  and  a 
metered  flow  of  the  additive  injected  with  an  electric  potential  of 
12,000  to  18,000  volts  between  the  injector  and  upper  honeycomb 
section.  Potentials  greater  than  18,000  volts  resulted  in  arcing 
between  the  fuel  injector  and  the  grid.  Atomization  benefits  of 
the  electrostatic  field  were  evident,  but  an  even  dispersion  across 
the  flame  holder  area  was  not  obtainable.  At  additive  flow- rates 
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greater  than  5#  by  weight,  flame  quenching  or  arcing  was  experienced. 
Elevation  of  the  inlet  air  temperature  to  500 *P  and  modifications 
to  the  injection  procedure  did  not  resolve  the  problem. 


Test  Results 

Flame  stability  limits  obtained  in  these  tests  are  presented  in 
Figure  4  and  Figure  5«  The  vertical  ordinate  in  both  figures  is 
fuel  equivalence  ratio  or  the  operating  fuel-air  ratio  divided  by 
the  stoichiometric  fuel-air  ratio.  In  Figure  4,  the  maximum  stable 
burning  limits  are  plotted  in  terms  of  burner  air  flow  velocity  at 
the  inlet  air  temperature.  Figure  5  presents  the  same  data  on  the 
basis  of  total  gas  flow  velocity.  While  the  velocity  limit  values 
differ  in  the  two  figures,  the  relative  ratings  of  the  additive 
effects  are  not  changed. 

Due  to  the  variations  in  inlet  air  temperature  found  necessary  for 
examination  of  these  additives,  the  influence  of  the  gaseous  addi¬ 
tives  upon  the  maximum  stanle  ourning  limits  is  more  easily  compared 
in  Figure  6.  Here  the  percentage  improvement  in  blow-out  velocity 
is  shown  for  each  addi tive  over  that  for  the  ammonia-air  flame  at 
the  same  air  inlet  temperature  and  fuel  equivalence  ratio.  The 
influence  of  pre-heating  the  ammonia-air  mixture  to  300*F  1  3  also 
compared  with  77*F  ammonia-air  limits  in  this  figure.  Also  plotted 
in  this  figure  is  the  improvement  in  flame  stability  limits  afforded 
by  the  additions  of  5#  by  volume  of  acetyler,.^. 
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The  following  observations  can  be  made  from  the  Figure  6  summarizing 
presentation.  In  5#  by  volume  additions  to  the  fuel,  none  of  the 
additives  tested  in  this  project  appear  as  beneficial  in  Improving 
the  ammonia  air  flame  propagation  rate  as  partial  pre-cracking  (28j<) 
of  the  ammonia.  Acetylene  additions  at  5#  volume  are  also  competitiv 
with  the  additives  investigated  herein.  The  flame  propagation  effect 
of  the  additives  varied  significantly  over  a  relatively  narrow  range 
of  fuel-air  ratio.  The  "fuel  type"  additive  (CO)  exhibited  a  marked 
improvement  in  the  fuel-lean  are?  only,  while  the  nitrogen  oxides 
effects  were  larger  in  the  fuel-rich  region.  Of  those  materials 
examined,  nitric  oxide  appears  most  promising  as  to  flame  propagatior 
rate  improvement. 

Quantitative  data  was  not  obtained  for  ammonium  nitrate  additions 
due  to  the  atomization  and  additive  distribution  difficulcies 
previously  described.  Qualitatively,  it  appeared  that  the  ammonium 
nitrate  might  be  improving  the  propagation  rate  of  the  flame.  When 
the  ammonia  flame  initially  separated  from  the  flame  holder,  a  pale 
green  reaction  zone,  characteristic  of  ammonium  nitrate  flames, 
remained  attached  to  the  flame  holder.  The  injection  operating 
difficulties  and  the  transient  character  of  this  condition  prevented 
its  detailed  examination. 


Burner  system  fouling  encountered  during  operation  with  NO  and  N0o 
is  also  of  significance  to  considerations  of  practical  application 
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of  these  materials  in  ammonia  fueled  engines.  In  our  tests,  the 
nitrogen  oxides  were  admitted  to  the  air  stream  prior  to  mixing 
with  the  fuel.  Downstream  of  the  air- fuel  mixing  region  is  a 
sintered  stainless  steel  filter  disc.  Deposit  formation  was  en¬ 
countered  where  the  fuel  and  additive  containing  air  flow  joined 
in  sufficient  quantity  to  clock  the  sintered  filter.  Chemical 
analysis  indicated  these  deposits  to  be  ammonium  nitrate.  Removal 
of  the  sintered  filter  screen  and  the  300°F  inlet  air  temperature 
used  for  N02  testing  alleviated  the  problem  in  our-  case,  but  this 
condition  could  be  detrimental  in  extended  engine  operations. 

While  only  comparative,  the  relative  values  reported  herein  on 
additive  effects  are  believed  to  provide  a  valid  screening  of 
gaseous  additives  based  upon  previous  correlations  between  this 
test  and  engine  performance.  However,  extensive  modifications 
to  the  bench  burner  approach  appear  to  be  necessary  before  meaningf 
data  can  be  obtained  on  materials  such  as  ammonium  nitrate  which 
cannot  be  readily  vaporized. 
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AMMONIA  FUELED  SPARK  IGNITION  ENGINE 
DISSOCIATOR  DEVELOPMENT 

Allison  Division  of  Generi 1  Motors  Corp. 

AMMONIA  D IS SOC IA T ION_BA CKGROUND  DATA_ 

Catalytic  Material  Investigation 

Early  in  1964,  Allison  made  a  survey  of  the  use  of  catalysts  in  the 
ammonia  industry.  As  a  result  of  this  survey,  the  following  four 
catalysts  were  selected  for  investigation: 

Triplv  promoted  iron  (ICI  35-4) 

Nickel  on  alumina 

Iron  on  silica  (G-47) 

Platinum  on  alumina 

Preliminary  teats  were  run  to  evaluate  the  effectiveness  of  the 
various  catalysts  vs.  no  catalyst  st  ammonia  flow  rate*  of  1.3  lb, /hr. 
and  at  various  temperatures.  Results  are  shown  ii«  Figure  1.  As  indicated, 
the  promoted  iron  was  the  only  catalytic  material  shoving  any  measurable 
increase  in  dissociation  rates.  No  further  tests  were  made  with  the 
other  three  catalytic  materials. 

In  August  1964,  a  subscale  ammonia  diosociator  was  designed  and 
fabricated  tc  determine  the  performance  of  a  dissociate'  jnder  simulated 
engine  operating  conditions.  Dissociation  of  the  am®  -  fuel  was 
achieved  by  passing  superheated  ammonia  vapor  through  a  bed  of  triply- 
-promoted  iron  oxide  catalyst  which  was  heated  by  engine  exhaust  gases. 

The  effect  of  exhaust  temperature  on  the  dissociation  level  at  various 
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ammonia  flow  rates  was  determined  by  operating  the  dissociator  at  exhaust 
gas  temperature*  above  and  below  the  equivalent  engine  temperatures. 
Results  of  this  investigation  are  shown  in  Figure  2. 

L»14I  Engine  Dissociator  System 

In  the  early  stages  of  the  present  contract,  and  before  actual  engine 
data  was  available,  Allison  was  given  the  assignment  to  design  end  develop 
a  dissociator  .for  the  L-141  engine.  Th<  following  design  criteria1,  /ere 
established: 

Fuel  enrichment  -  2.5  percent  Hydrogen  by  weight 
XH2  ■  H,  x  100 

H24N2+  NH3 

Fuel  temperature  -  The  partially  dissociated  ammonia  fuel  to  be 
cooled  to  a  temperature  compatible  with 
requirements  of  tha  carburetor. 

Physical  sire  -  To  permit  installation  in  the  M-151,  1/4  Ton  vehicle 

Design  parameters  -  Based  on  2400  r.p.m.  gasoline  engine  deta 

(See  Teble  I). 

Dissociation  Producto  Cooler 

The  dissociation  products  cooler  -  ammonia  vaporiser  consisted  of 
two  helical  coils  enclosed  in  a  3-1/2  inch  stainless  steel  shell.  The 
high  temperature  partially  dissociated  ammonia  from  the  dissociator 
was  cooled  as  it  passed  through  the  inside  of  the  coiled  tube  and  the  heat 
rejected  was  used  to  vaporize  the  incoming  liquid  ammonia  which  passed 
over  the  outside  of  the  tubing.  Calculations  based  on  tests  of  the  sub- 
seal*.  dissociator  indicated  that  the  dissociation  products  would  be  cooled 
to  approximately  80°  F,  The  heat  rejected  by  cooling  tha  dissociation 


Subscale  Dissociator,  Dissociation  Versus  Exhaust  Temperature 
0.72U  Pound  ICI35-U  Catalyst. 
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TABLE  I 

Dissociator  Design  Parameters  at  Selected  L-141  Engine  Condition 


Engine  speed 

Manifold  pressure 

Ammonia  flow  rate 

Air  flow  rate 

Exhaust  gas  temperature 

Ammonia  inlet  temperature 

Ammonia  pressure 

Dissociator  products  temperature 


2400  r.p.m. 

18"  hg.  a.b.s, 
40.35  lb. /hr. 
245  lb. /hr. 
1275°  F. 

32°  F. 

1  58.!?  p. 8.1.4. 
o 

cooler  inlet  820  p. 
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products  from  820°r'  to  80°F.  was  sufficient  to  vaporize  about  b3?«  of 
the  incoming  liquid  ammonia.  The  overall  size  of  the  dissociation 
products  cooler-ammonia  vaporizer  was  3^-inch  diameter  by  "J.  inches 
in  length.  Design  details  are  shown  in  Figure  3,  Allison  drawing 
EX-73993*  Continental  drawing  595713, 

Preheater-Pi  3socit\tor 

The  dissociator  unit  for  the  L-llil  engine  was  designed  to  be  attached 
to  the  engine  exhaust  manifold  and  to  utilize  exhaust  waste  hed.  to 
preheat  and  dissociate  the  ammonia.  The  preheater  section  consisted  of 
two  helicol  coils  of  £-inch  stainloss  steel  tubing  wound  in  opposite 
directions  and  joined  by  a  reverse  band  to  minimize  thermal  distortion, 
catalyst  bed,  fabricated  of  2  3/8-inch  stainless  steel  tube,  was 
positioned  inside  the  preheater  coil.  Eight  longitudinal  fins  on  the 
catalyst  and  supporting  the  preheater  coil.  Five  pounds  of  triply-pro¬ 
moted  iron  oxide  catalyst  (ICI  35-lt)  were  packed  iu  the  catalyst  tube. 
The  catalyst  was  retained  by  cones  of  20-mesh  stainless  steel  wire 
cloth.  The  preheater  coil  and  catalyst  tube  were  encased  in  a  Uj-inch 
steel  shell,  1?  inches  long.  Diametrically  opposed  exhaust  ports  were 
located  at  opposite  ends  of  the  shell.  The  exhaust  gas  flowed  axially 
through  the  shell,  acrors  the  preheater  coils  and  parallel  to  the 
catalyst  tubs. 

The  dissociator  was  designed  to  complete  vaporisation  of  the  ammonia, to 
hast  the  ammonia  vapor  to  dissociation  temperature,  and  to  dissociate 
1U,2  percent  of  the  ammonia  (providing  2,5  percent  hydrogen  by  weight). 
During  the  single  pass  through  the  preheater-dissociator,  the  exhaust 

gas  was  cooled  from  1275°  to  975°?.  Design  details  of  the  preheater- 
dissociator  are  shown  in  Figure  it,  Allison  drawing  SK-15363,  Continental 
drawing  595711, 


gap  ttepo: 
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Two  dissociator  systems  were  fabricated;  one  for  bench  testing  at 
Allison,  and  one  for  engine  testing  at  Continental. 

Preliminary  Testing  of  L-141  Dissociator 

The  L-141  dissociator  system  consisting  of  the  EX-73993  cooler- 
vaporizer  and  SK-15363  preheeter-dissociator,  was  installed  in  the 
Allison  burner  teat  rig.  Liquid  ammonia  was  available  in  the  test  cell 
at  room  temperature  and  at  225  p.s.i.g.  pressure.  Hot  exhaust  gas  was 
furnished  by  a  JP-4  (hydrocarbon)  fueled  combustor.  Instrumentation  was 
provided  to  measure  presares,  temperatures  and  flows  at  all  critical 
locations. 


Performance  data  obtained  at  temperature  and  flow  conditions 
corresponding  to  engine  operation  at  speeds  from  800  to  3600  r.p>m.  and 
at  manifold  pressures  from  12  to  30  inches  of  mercury  are  shown  in 
Table  II,  The  fuel  enrishment  provided  by  the  dissociator  at  various 
simulated  engine  conditions  is  shown  in  Figure  5. 


The  excessively  high  cooler  outlet  temperatures  shown  in  Table  II 
were  due  to  liquid  ammonia  being  thrown  to  the  outside  of  the  coiled 
passages  and  a  large  amount  of  the  heat  transfer  surface  was  not  being 
effectively  used.  Th*8  situation  was  corrected  by  replumbing  the 
cooler  so  that  the  liquid  ononis  was  inside  the  tubes  end  the  hot 


produets  were  on  the  shell  side. 


1S&S-2 

IgULBisaoelator  Svoteg  Test  Data 
SlmtXAted  Bnglne 
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Preliminary  engine  resting  of  the  dissociator  system  was  terminated 
at  an  early  date  because  of  Id  35-4  catalyst  disintegrated  into  a  fine 
powder  that  could  not  be  retained  in  the  catalyst  bed. 

Modified  L-141  Pissociator  System 

The  tendency  of  the  granules  of  ICI  35-4  triply  promoted  iron  catalyst 
to  disintegrate  under  temperature  and  vibration  conditions  encountered 
in  engine  operation  rendered  it  clearly  unsuitable.  Several  attempts  were 
made  to  use  this  catalyst  in  another  form,  such  as  sintering  into  a 
solid  plug  and  plasma  spraying  on  stainless  steel  mesh.  These 
approaches  were  all  unsatisfactory. 

A  number  of  other  catalysts  were  investigated  as  substitutes  for 
ICI  35-4.  Among  these  were  promoted  iron  (KM1R) ,  stainless  stell, 

Incomel  and  nickel.  PorouB  nickel  in  the  form  of  1  Foametel”  was 
finally  selected  due  to  the  large  active  surface  area  and  high  mechanical 
strength.  The  relatively  low  activity  of  the  nickel  catalyst  necessitated 
a  complete  redesign  of  the  prehe^ter-dissociator  unit.  A  new  set  of 
design  criterian,  based  on  actual  engine  test  results,  were  established 
as  follows: 

Fuel  enrichment  1.25  percent  hydrogen  by  weight 

at  1000  r.p.m.  full  throttle  to 
3.15  percent  hydrogen  by  weight 
at  4000  r.p.m.  full  throttle. 

Fuel  temperature  The  partially  dissociated  ammonia 

fuel  to  be  delivered  to  the  carburetor 
at  a  temperature  of  50  -  100  F. 
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Design  parameters 


Based  on  2400  r.p.rn. ,  full  throttli 
operation  of  the  L-141  engine  on 
hydrogen  enriched  ammonia  fuel 
(see  Table  III). 


TABLE  III 


Revised  Design  Parameters  for  L-141  Engine  Dissociator 


Er  gine’  Speed 
Manifold  pressure 


Ammonia  flow  rate 


Air  flow  rate 


Exhaust  gas  temperature 
Ammonia  inlet  temperature 
Ammonia  pressure 


2400  r.p.m. 
29"  hg  a.b.s. 
37  lb; /hr. 

220  lb.  /hr. 


1235°  F. 
32°  F. 


53. 5  p.Sci.a. 


Dissociation  products  temperature  at  catalyst  outlet  1025°  F. 


The  preheater-dlssociator  was  redesigned  with  an  automatic  bypass 
arrangement  so  that  only  part  of  the  ammonia  vapor  would  flow  over  the 
catalyst  bed.  In  order  to  obtain  a  low  pressure  drop  across  the  catalyst, 
the  Foametel  nickel  was  machined  into  16  discs,  each  29/32  inch  thicx  and 
3.275  inches  in  dismeter.  The  discs  were  separated  in  the  catalyst  tube 
by  3/32  inch  spacers,  and  internal  manifolding  was  provide!  so  that  ammonia 
flowed  through  the  16  discs  in  parallel.  The  catalyst  bed  was  installed 
in  a  stainless  steAl  tube  16”I/2  inches  long,  having  32  longitudinal  fins 
to  increase  heat  transfer.  Two  double  helicol  coils  of  5/16  inch* 
stainless  steel  tubing  were  wound  concentric  with  the  catalyst  tube, 
providing  four  parallel  flow  paths  for  the  aarmoniu  vapor.  A  thin 


cylindrical  shield  was  placed  around  the  catalyst  tuba  to  dir  act  the 


CAE  Report  No.  105U 
Appendix  III 
Volume  II 


-lli- 


exhaust  gas  flow.  The  preheater  coils  and  catalyst  tube  were  encased  in 
a  6  inch  stainless  steel  shell,  a  split  3-1/2  inch  tube,  extending 
the  length  of  the  preheater-dissociator  shell ,  contained  a  70  inch  length 
of  1/2  inch  tubing  which  served  as  an  auxiliary  evaporator. 

Exhaust  gas  front  the  engine  enters  through  the  center  of  one  end 
of  the.  preheater-dissociator  shell,  flows  along  the  catalyst  tube,  back 
across  the  preheater  coils,  out  the  side  of  the  chell  into  the  auxiliary 
vaporizer  section,  axially  through  the  auxiliary  vaporizer  and  exits  into 
the  vehicle  exhaust  system.  Design  details  of  the  modified  preheater- 
dissociator  are  shown  in  Figure  6,  Allison  drawing  EX-79016,  CAE 
drawing  595712. 

Electrically  Heated  Auxiliary  Dissociatcr 

When  the  L-141  engine  is  operating  under  part  throttle  conditions  with 
ammonia  -apor  alone  as  the  fuel,  the  exhaust  temperatures  are  below 
1000°  F.  Under  these  conditions,  the  nickel  catalyst  is  relatively 
ineffective.  It  was,  therefore,  decided  to  design  and  fabricate  an 
electrically  heated  auxiliary  dissociator  for  use  during  starting  and 
low  load  operation. 

The  size  and  geometry  of  the  auxiliary  aisaociator  was  determined 
primarily  by  available  heater  and  tubing  sizes.  The  auxiliary  dissociator 
consists  of  a  concentric  tube  heat  exchanger  with  a  Foametai  nickel 
cylindrical  catalyst  and  an  electric  cartridge  heating  element.  A  5/8 
inch  diameter  by  6  inch  long,  28  volt,  450  watt  Firerod  cartridge, 
manufactured  by  Watlov  Electric  Manufacturing  Compauy,  was  selected  for  use. 
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EX79016  Preheater-Dissociator, 
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The  Porous  nickel  catalyst  cylinder  which  surrounded  t-'s  heating  element 
van  1,5  inches  in  diameter  by  5/16  inch  wall  bj  7  inches  long  and  weighed 
approximately  1  pound.  This  was  encased  in  a  2 '1/4  inch  tube  by  6  inches 
long. 


Preheated  ammonia  vapor  from  the  dissociator  outlet  enters  the  outer 
annulus  of  the  auxiliary  dissociator.  Approximately  2/3  of  the  rntering 
amrtMnia  is  bypassed  directly  to  the  auxiliary  dissociator  outlet  through 
a  hole  drilled  in  the  inner  tube.  The  remaining  ammonia  flws  axially 
through  the  t  .er  annulus,  back  along  the  heating  element,  then  radially 
Chrougn  the  catalyst  cylinder.  From  the  auxiliary  dissociator,  the  disso¬ 
ciation  products  are  routed  into  the  cooler  side  of  the  EX-73993  vaporizor 
cooler.  A  cross-sectional  sketch  of  the  auxiliary  dissociator  is  shown 
in  Figure  7. 

JJiE'iSE  Rig  Testing  of  the  Redesigned  L-141  Dissociator  System 

The  redesigned  L-141  dissociator  system,  consisting  of  the- 
EX-73993  cooler-vaporizor,  the  EX-79016  redesigned  preheater-dissociator 
and  the  electrically  heated  auxiliary  dissociator,  together  with  the 
necessary  plumbing  and  instrumentation,  was  installed  in  the  burner  test  rig,  A 
schematic  of  the  test  setup  is  shown  in  Figure  8.  Performance  data 
obtained  at  temperature  and  flow  conditions  corresponding  to  engine 
operation  at  speeds  from  200  to  4000  r.p.m.  and  at  manifold  pressures 
of  15  and  29  inches  of  mercury  are  shown  in  Table  IV.  The  fuel  enrich¬ 
ment  provided  by  the  dissociator  at  various  Simula tad  engine  conditions 
with  the  auxiliary  dissociator  off  is  compared  to  the  desired  and 
predicted  enrichment  in  Figure  9. 
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Fig*  Burner  Test  Rig  Setup  for  Redesigned  L-li*l  Disscciator, 
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Engine  testing  of  the  redesigned  L-lhl  dissociator  system  is 
discussed  in  the  main  body  of  the  report* 

Design  calculations  for  the  original  L-llfL  engine  dissociator 
system,  the  redesigned  system  using  nickel  catalyst,  and  the  electri¬ 
cally  heated  auxiliary  dissociator  are  included  as  Appendices  kt  B  and 
C. 
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Calculated  air  and  fuel  flows  for  L-141  engine 
Engine  displacement  —  141  ln^ 

Swept  volume  per  revolution  -  g  ~  *0408  ft^/rov 

Density  of  airi 

(assume  T  “  110°?  -  570°H) 


AJ2S* 


,  00232  /° 


Density  of  ail*  at  various  manifold  pressure* » 


/° 

S'/'/? 


*.00232X/2  * 

*00*39  X  /&  - 

*00232X34  ~ 

<00232X30  - 


£6/&S 
.0V'?  H 
t0S27  " 

;  0696  " 


Air  weight  flow  per  revolution  at  various  manifold  pressures, 
volumetric  efficiency)* 

e* 

(Assume  lOOJl 

.6Y06X.0273 

*  ,007/4 

Ib/rev 

'<?*  Jt? 

•0400X07/7 

-  .00/70- 

lb/rev 

J'/'S 

.0'/4?x.o%rr 

-  ,60727 

lb/r®* 

fo'//? 

«  .00234 

lb/rev 

Assume  air-fuel  ratio  6,07i1. 
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ondllions  Selected  for  DesignJ^oink 


Ritt 

Munifold  pressure 
Fu*l  flow 

Air  flow 

Exhaust  mass  flow 
Exhaust  gas  temperature 
Ammonia  temperature 


Ammonia  pressure 


Dissociation  products  temperature 
at  cooler  inlet 
(for  2.5S  i\ydrogen  enrichment 
and  5  lbs  ICI  35-4  catalyst) 


2400 
18”  hg 

40.35  lbs/hr 
245  lbs/hr 

285.35  lbs/hr 
1275°F 

32°F 

58.5  P*i* 

820°F 
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MH3  flow _  _  lbe/hr 

Catalyst  weight  lb 
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Configuration:  Double  helix  heat  exchanger 
Liquid  MH3  on  outside  of  coil 
Dissociation  products  inside  coil 


Assume  18  turns  £  OD  x  .032  wall  tubing  each  helix  length  of  tubing 
(2.76  7J  x  18)  +  (1.6  ttx  18}  =  24S  in 
Surface  area  . 5  ff  x  248  =  39C  in2  =  2.7  ft2 


Cold  side  conditions:  MK3  flow  rate 

40-35  lbs/hr 

Pressure 

58.5  psia 

Inlet  temperature 

32°F 

Outlet  temperature 

32°F 

Hot  side  conditions:  Gas  flow  rate 

40.35  lba/hr 

Pressure 

58.5  psia 

Inlot  temperature 

820°F 

Cooler  -  hot  side: 
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~  39, 006  44>/drjP2 

^  _  Mo/  &  _  >6343  X  39000  m 


„  MMJ4L  44/Mr 
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For  exhaust  gas  flow  over  outside  of  tubes 


Flow  area  =  cross  section  of  outer  shell  (inside)  -  cross  section  of  catalyst 
tube  (outside)  -  area  of  preheater  tube  (end  view) 
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Dlssociator  Section  Heat  Balance . 
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Design  Calculationa  for  Radeaignad  M41  Diiaociator 
(EX-79016  with  Nickel  Catalyet) 


iOXii&SZX.  e°f  tjctton 

(QJ  jf  30  /,0L?S  6LVl/f  /  70  m  Jzxp  Jl?  . — . 

_ Jpa/. -  f/ry.70.-.  //0m*. 70X731. 

~e& fx£ 7s'a^Jtoa3&0XI . 7 '"777777 7 777 


XJ2Z 


/5?Za/~ _ 2SL  ^ 


Preheater  Section 

/f  00  X.  073  WO/T  ff  7#0u*tp . .  .  • . 

..  70  &rMS.  3.7373  3>0  +  _  70  &J073 /&  . -  . 

_ . ps&r.xyt^ 

. x#iz/  3xm/o0<?-.  3/JSjrX7370.  z/pj&m*  *-<£X3y2  7  . 


Catalyat  Section 

*\  ........  ..  . 

.  3,77  00/  .073  waSS  X/Xjf  m/4  .37  j&*S. 

......  009  X  .070  /  /3f  .  .  .  - - - -  --  - 

Xu3p  su*y&c3  .  „3*77/r.X  /&S1*  / 7 77*  23..X 03. // * - 

33X3/  .79  X  /S33^3$3^73  703/*.* 
S^&J' JUtryilC*  733+  702  - .... 3.30 X7^..  ... .  ...  - . 


CAS  Reoort  Ho.  105k 
Aopendix  III 
Volume  II 


Catalyst 


{<££)..  VS%  . .  . 

//  aiscs ,  J.J7S  x  .9a£  Jfcat  . . 

.  _ _ _  _ _ _ ....  . . .. 

of di&/ys£ . . .  ._  _ _  ■ 

//tttb3x  ^22/Z/uj3/  .VS  r  J/6.2  /X.  . -  . 


.....Omss  J&cihr??.  Xa&i&sX ^(pCiTeaC /a . ytfVs 

•  "Z£j±”*  ,  „  -  •  .  -* - 

-  *?06u^t  ._ J :. — f= — .S9 - 


l 


CJtE  Report  No.  105U 
Aopendix  III 
VoluinR  II 


-1*1- 


CalcuUted  Pre««ur«  Drop  in  Catalyat  ?oaa«t»l  Ntckal,  43X  Dan»ltff  29w  ha. 
4000  RPM  Engine  Conditions 


s O&qo  Tf&ztjDen f/bryj)  .... 


.//  ~ jQr&tJur? _ *0 josia,. - ^ .. _ 

_/?  . . . .  -  _  - 

.  ..  ms**  .  _ m.07sfj#  .  . 

.$ .  sJ%/X?rfaa/, /&s£<ufy.  _ *.  MJJ7T2/J&/7* 

J{  B/mcauTy.  .  .  . . . . */f//0*a//2jec.. 

$c  =-  tifctot&t&ataft coiisfarU _ &SJ./7 j#/sacL  _ 

M*.  /friol&Cu&r  ut&p££. ... _ _  _ _ +S7JZ/&4M& 

7?  -  <?as  C&ulT&hJ. . . . */W6J7&/&m£ 

.  T  -  ^yo<?ra7Z>r<?  .  /Wf"  s  /704  V? . 


.; /fam  $7&9in  f.  Zlcav&z  ( _ _ 

.  a*  ~-2X/aetot~2 X  /vy. _ 

./? •  .  /'SMWj  /?$ . » . AfaX/0*Jt~*  _ 


CAE  Reoort  No,  105U 
Appendix  III 
Volume  II 


-i»2- 


J 


-  texwf-  /$*  _  py.tiry/o*-/*!* 

X  "  .07SS-  “  .07?r 


u 

y 

ii 


Xa/pr^tG  2t2&8X/aJ£>//SMX/7MX/S''sC*X.MS‘27 

Mfr  *  /7 X  32.//  ~~  x  26.3  x  /C* 

/?*  Zfo  ft  '  /.<S>2/76c-X:o7srl*j^  A  #6  X /o  ' 


2/P7&  9  2 X  7SV&_  X  7702  Xf.MS27? 

Mg?  ~  77  X  22.77  ‘  .^7r 

Su&sfxfu/esn# : 

7Y.62X70*  -/$*  <  /  ,  , 

.&7S2  ~  7X2X70  +  (/£0X/0^  /  .  2670") 


/§  -  7yJ?X/6(  -  2.63X70-6  s  72.4?  X7Cc 

/?  =  f6//f2  =  jT^  yO^Cc 


Pxessi^r  dPrzgo  ~ 


.  ! 
.  / 


I 

J 

if 

] 

I 

i 

I 


k 


CAE  Report  No.  10£i* 
Appendix  III 
Volume  II 


Calculated  Praaaura  Prop  In  Preheater  Coil  29"  hg,  4000  RIM  Inline  Ccndltlona 

of  n&a&J  »3/0S- .070  *.2¥2Su*t  *  .0042  ft 

Jh?<?  zrsJec&w  feer  *_ .  00033  /&  * — 


@  =  7X00033  ft*  *  -  te'Soo  &/&-#*■ 


A/  ?*L£  -  l  X/3000& 

A  “  .  OSS  06/A/& 


't - 


7620 


/£  p  ££/£  *  A  -  - - 

-Ulf  -  y  *  *  ^  ^  -  .A7//7  &/S0C. 

/&  _  V?A/7*01 f  --  ■  ■  ~  -  •  . 

-  00/%4*>*x/yy  _/7.y/7y& 

::.j;:i::_  j_r_z_ 

Jfer  j&r&L M&n? .afrqto  tf&u  ol  yo^oe ...  (ftonfa}. _ 

4^  “  /7Y'2.  \f  or/  tC/o/'**  . . . .  . 

.  (TrforAsJ  _ 


4A - .  ms x . m x(.00//rfxJ7.X x  M- /  Xfim*)*- 

.  -  /r?‘tx.ayi  /  /j.jsx/fjr/xyx  ~ 

. =  .  ..Jl5L-/>£i _ - _ 


CAE  Report  No,  10$U 
Appendix  III 
Volume  II 


■feU- 


StresG  in  Catalyst  Tube  at  60  pals 
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Appendix  III-C 

Design  Calculations  for  Electrically  Heated  Auxiliary  Dissoclator 

Design  Point:  1200  RPM 

15"  hg.  Manifold  Pressure 

810°F  Exhaust  Temperature  (Fig.  (M) 

10  lb/hr  NH3  Flow  (Fig.  C-2) 

Assume  desired  enrichment  =  17.  H2  by  weight 

r  5.657,  NH3  dissociated 

Assume  NH-j  vaporized  and  preheated  to  700°F  in  existing  dissoclator. 

(Figure  C-3) 

Fron^  catalyst  performance  curve  'Figure  0-5) ,  167  dissociation  may  ba  obtained  at 
a  catalyst  bed  temperature  pf,  1000°F  with  a  flow  of  1.4  Ibs/hr  through  .05  lb  nickel 
on  alumina  catalyst. 

Required  flow  thru  dissociator  -  35.27. 

or  .352  x  i0  lb/hr  =  3.52  Ib/hr 

required  to  heat  35.27.  of  Wf  from  700°F  to  1000°F 
Qh  »  3.52  lb/hr  x  .61  8tu/ibl,F  x  300°F  *  644  Btu/hr 

required  to  dissociate  5.657  of  Wf 
Qj  =  1354  Btu/lb  x  10  Ib/hr  x  .0565  =  765  Btu/hr 

Total  Q  =  644  «■  765  =  1409  Btu/hr 

1409  Btu/hr  x  .2931  watts/Btu/hr  =  413  watts 

Assuan*  24  volt  source, 

24  volts 


17.2  amps 
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it  700°F  Inlet  temperature,  csaua*1  15  psia  pressure 

.63T  .63  x  1160  ,  „  „ 

v  = - *  =  - - — - —  -  48 v8  ft- /lb 

P  15 

C  3<>52  lb/hr  flow  thru  ditsociator 

V  =  3.52  lb/hr  x  48.8  ft3/lb  *  172  ft3, 'hr 

*  2.9  ft3/win 

rom  catalyst  performance  curve  (Figure  3)  dissociator  with  nickel  on 
lualna  catalyst  will  be  1.5  in.  long  x  1.25  in.  diameter  with  .05  lb. 
italyst,  1.4  lb/hr  NH3  flow. 

seizing  for  3.52  lb/hr  flow: 

=  2.5  x  size 


ea  of  1.25  in.  diameter  catalyst  bed  =  1.227  in.' 


1.227  in.2  x  2.5  3  3.06  in.2 


2"  diameter 


auired  catalyst  bed  6ize  1.5  in.  long  x  2  in.  diameter 
talyst  weight  *  2.5  x  .05  lb.  =  .125  lbs. 


rting  conditions  (at  59°F  Ambient  temperature) 
200  RFM  cranking  speed 
29"  hg.  manifold  pressure  (14.2  psia) 


ume  NH3  vapor  temperature 
aity  of  air 


59°F  (to  carb.) 


e  -  f= 

Ki 


14.2  psia  x  144  in.2/ft2 

53.3  it. lb./lb°R  x  519°R 


«■?/.  su/<*3 

e  v  *  -»  aw:  a  v 


sity  of 


14.2  psia  x  144  io.?/ft2 
90.77  ft. lb/ib°R  r  519°R 


.0435  lb/ ft3 
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X  H2  by  weight  =  1_  x  .072  W  .012  *  1.2X 

Cc'd  start  condition  (Assuan*  -40°F  Ambient  Temperature) 

200  RPM  cranking  speed 
29"  hg.  manifold  pressure 
4.48  lb/nr  NH3  Flow 

Qh  required  to  raise  temperature  of  4.48  ibs/hr  liquid  from  «40°F  to  -3°F 

Qh  =  4.48  Ib/hr  x  1.065  Btu/ib°F  x  35°F  «  167  Btu/hr 

Qv  required  to  vaporize  4.48  Ib/hr 

Qv  =  4.48  Ib/hr  x  589.3  Btu/lb  *  2640  Btu/hr 

Total  Q  *  2640  +  J 67  =  2807  Btu/hr 


Assume  regenerative  system  to  cool  NH3  out  of  dlssociator  from  1000°F  to  100°F 
Q  available  from  regenerator 


1.59  Ib/hr  x  .61  Btu/lb°F  x  900°F  «  872  Btu/hr 


Net  Q  requirement  2807  -  872  =  1935  Btu/hr 

1935  Btu/hr  x  .2931  watte/Btu/hr  *  566  watts 


566  watts 
24  volts 


23.6  amps 


Total  powe^r  requirement  during  cold  start  (*40°F)  vaporiser  plus  dlssociator 
413  ♦  533  *  946  wabts 


17.2  ♦  22.2 


39.4  amps  (at  24V) 


Valght  of  24V  Battery  with  capacity  of  54  amps  for  5  minutes  *  34  lbs 
(Exide  #53033) 


•73- 
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Check  Point:  800  RPH 

22“  hg.  Manifold  Pressure 
675°F  Exhaust  Temperature 
7.5  ib/hr  N43  Flow 


Flow  thru  dtsaociator  =  .352  x  7.5  Ib/hr  *  2. 64  Ib/hr 

From  figure  C-3,  NH3  is  preheated  to  50()®F 
Q  available  1409  Btu/hr 


Assume  (for  trinl)  5.65Z  dissociation  of  7.5  Ib/hr 
(or  16Z  dissociation  of  2.64  ib/hr) 

'or  155  H2  by  weight) 

Qd  =  1154  Btu/lb  x  7.5  ib/hr  x  .0565  =  573  Btu/hr 

leaving  1409  -  573  -  836  Btu/hr  for  heating  NH3 

Temperature  in  catalyst 


500 


#  836  Btu/hr _ 

2.64  lb/hr  x  .61  Btu/lb°F 


1020°F 


which  is  the  correct  temperature  for  this  dissociation  level  'Figure  C-5) 
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Percgnt  Dla«oclatad  Veram  Agiaonlt  Teoperature 
fet ' fcickel  on  Alumina  Catalyst  *  Reference  1 


800  900  1000  1100  1200  1300 


Ammonia  temperature  --■°F 
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ABSTRACT 

A  stu.ty  was  conduct**!  to  «le*t •  •  r'n f n>'  the*  feasibility  of  u«tng  two 
known  sources  of  high  energy  to  overcome  some  of  the  dlff lenities  In¬ 
volved  tn  using  ammonia  ns  a  reciprocating-engine  fuel.  The  two  approaches 
were:  (a)  The  employment  of  a  high  voltage,  high-energy,  spark  disehirge 
across  as  large  a  gap  as  possible  In  the  combustion  chamber.  "Bomb" 
tests  indicated  that  the  desirable  combustion  characterlst les  could  be 
realized  and  that  combustion  of  the  mixture  was  reliable  under  conditions 
of  temperature  and  pressure  which  would  cause  a  standard  ignition  system 
to  fail  completely,  (o)  The  dissociation  of  inlet  air  (nitrogen)  by  a 
strong  local  radio- frequency  field,  adding  the  recombination  energy  to 
the  cyl inder-t rapped  gases.  Differential  equations  of  reaction  rates 
were  solved  b>  using  computer  subroutine  PRECOR.  Product  lifetimes  were 
less  than  desired,  but  temperature  dependent. 


KEY  WORDS 


Fuel 

Energy  Depot 

I  nr i f Ion 

Ammonia 

Dissociation 

Engines  (Reciprocating) 


v 


R- 1122 


CAE  Report  No.  1054 
Appendix  IV 
Volume  II 


TABLE  OF  CONTENTS 

ABSTRACT . v 

BACKGROUND  .  . .  I 

HIGH-ENERGY  DISCHARGE  STUDY  .  2 

Apparatus  . 3 

Test  Proftram . 4 

Resu'ts  .  6 

Cunt  :a*>  ions  .  . .  8 

Recommendation .  r 

RAOIO-FRECUCNCY  DISSOCIATION  STUI-Y .  10 

Kinetics  of  the  Qlssociated-Mi trogen  Air  System  ........  12 

Primary  Reactions  .  12 

Secondary  Reaction; . 13 

Other  Reactions  . .  17 

Reaction  Rates  .  18 

Table  I,  Reaction  Rates  tor  the  Oi  ■  sot  iutod-Ni  t  rogrn 

Air  System  (Pri-wiy  Reactions)  .  K’ 

Table  2,  Reaction  Rates  for  the  Oi ^mx- iot ed-N i roqen 

Air  System  (Svcor.laty  Reaction*)  . .  20 

Tyne  a .  20 

Type  b . 21 

Type  t .  22 

Type  <1 . 23 

Table  3,  Unimportant  Reactions  for  the 

Oissocialevt-IJitro*t*  i-  oil  Systi  n .  24 

Cu  liu  ie  l  ion  of  Nitrogen  Recoil,  i oat  ion  Times  .  25 

Mode  i  I .  25 

Model  2 .  26 

Model  3 .  28 

Model  4 .  30 

Cone  In  - ions  and  Recommrndat ions  .  33 

REFERENCES .  35 

FIGURES  (list  on  p.  vill .  38 

APPENDIX  (Subroutine  RREC0R) . . .  A- 1 


R- 1 122 


CAE  Report  No.  1054 
Appendix  IV 
Volume  II 


Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Fig. 

Pig. 

Fig. 

Pig. 

Pig. 

Fig. 

Fig. 

Fig. 


LIST  OF  FIGURES 

1  Diagram  of  Test  Setup,  High  Energy  Discharge 

2  High-Energy  Pufser  Discharging  in  Air 

3  Two  Electrode  Configurations,  Showing  the  Tendency  to  Arc 

to  the  Cylinder  Head 

4  Sample  Data  From  Oscilloscope  Camera 

5  Chamber  Assembly  Drawing 

6  Top  Plate,  Type  1 

7  Top  Plate,  Type  2 

8  Combustion  Time  Versus  Pressure  Percent  of  Hydrogen 

9  Combustion  Time  Versus  Total  Pressure  for  Hydrogen 

10  Combustion  Time  Versus  Discharge  Energy  for  Hydrogen 

11  Combustion  Time  Versus  Pressure  Percent  Ammonia 

12  Combustion  Time  Versus  Total  Pressure  for  Ammonia 

13  Combustion  Time  Versus  Discharge  Energy  for  Ammonia 

14  Combustion  Time  Versus  Temperature  for  Ammonia 


15 

Pressure 

Rate  Versus  Temperature  for  Ammonia 

16 

Model 

1. 

10%  Initial 

% 

Dissoc iat ion 

17 

Model 

2, 

10%  initial 

No 

Dissociat ion 

18 

Mode] 

3, 

1%  Initial 

N, 

Dissociation 

19 

Model 

3. 

5%  Initial 

"a 

Dissociat ion 

20 

Model 

3, 

10%  Initial 

N 

Dissociation 

21 

Model 

3, 

1%  Initial 

Na 

Dissociation  600°K 

22 

Models  3 

9r.d  4,  10% 

Initial  Na  Dissociation 

A- 1 122 


CAE  Report  No.  10r* 
Appendix  IV 
Volume  II 


BACKGROUND  NOT  REPRODUCIBLE 

As  <s  r.;  ii U  of  th?  Erwrjy  D«*pot  C-mr-p* ,  catting  for  in-ihe-f lel<J 
production  if  liUor*i«»I”C*'v^!>u-.t  i  m-'.-n.jH?.  ft*  the  feasibility  >if  u'iny 
anhy  Irnt*  a  n'aio  ar-  a  fuel  has  bornm*  tJf  Hi  *r«*st. 

The  suitability  of  nvioiiia  as  a  fuel  for  the  internal  combustion 
engine  depends  on  the  answers  .-•>  a  number  of  questions,  So»ie  of  these 
i  tvolve  the  co-»t  and  l*usi>i*ity  of  projucinq  ammonia  in  the  field; 
olhi-i  s  r*lnte  to  the  ciiaract-  i  i  .tics  of  an  a* t  >-nia-fu..*led  engin  ,  such 
as  o  a  <■  e  of  starting,  power  out  *ut ,  (ability  *f  operation,  e>.c.;  still 
others  concern  the  m  ><Ji  f  icat  ion-  that  might  bn  rouu*red  to  enable  the 
c fwjine  to  use  ammorf.i  as  a  fuel,  enoino  corrosion  resulting  from  the 
combustion  products,  and  so  or. 

Two  peculiarities  of  the  ciutinn  of  thi  >  gas,  as  compute*!  with 
c*«ni»u  .1  ion  of  mrnnl  hy  'roiarh-n  fuels,  are  of  particular  Importance. 

The  fir-*,  is  the  relatively  lar.jn  ...r  unt  of  inerjy  necessary  to  lynite 
the  I'uctarts,  an<:  the  second  i.  the  relatively  slow  speed  of  flame  prop- 
agat ion. 

It  was  therefore  proposed  to  dissociate  or  ionize  cither  the  in- 
tak*  gases  or  the  cy • indcr- trapped  gases,  to  aid  and  accelerate  the 
oxidation  pr«>cess  in  the  co>ibu-  lion  chamber;  and  this  pr«>gram  was  under¬ 
take  to  examine  the  feasibility  of  using  two  known  high-energy  sources 
to  achieve  these  ends  without  major  modifications  to  present  engines. 

Methods  similar  to  those  selected  had  already  been  used  on  a  larger 
scale,  to  produce  hi«,h  energy  for  mobile  radar  equipment  and  bigh- 
tempi*rdturo  fu;ing  eeuipm^rt. 

The  present  work  was  divided  into  two  basic  tasks:  (a)  the  inves¬ 
tigation  of  an  arc  discharge  of  muchhighcr  energy  than  is  custo*iary  in 
spark- ign i t ion  applications,  and  (b)  the  dissociation  of  molecular  nitro¬ 
gen  by  passing  the  intake  air  through  an  intense  local  radlo-freruicnry 
(RF)  field. 
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HIGH-ENERGY  DISCHARGE  STUDY 

It  has  been  shown  In  pr  'vl  u1  that  th<*  combu  M  *r*  f  atunii 

ca.'  b*.  initiated  thi--ujh  spark  devices.  Hn/ever,  twit  nary  •  park  devices 
operating  on  pure  ammmia  and  air  Mixtures  have  I  *ft  *»uch  to  b«s  dt?slred, 
since  the  oner  jy  they  Make  available  Is  marginal  for  Ignition.  When 
supplied  at  a  point  source,  the  energy  does  not  rnmper-ate  for  the 
extraordinari  ly  low  flame  speed  of  ammonia  fuel. 

Th<*  purpose  of  this  stu.iy  >«a$  to  learn  whether  or  not  a  high-energy 
line  discharge  by  the  longest  possible  path  could  alleviate  those  problems. 

It  was  first  assumed  that,  if  many  times  he  minimum  ignition 
energy  were  aval  iablc,  amnnia  decomposi don  would  (<»n  a  sufficiently 
large  scale)  result  in  the  realization  of  better  burning  propert ios. 
Secondly,  It  was  assumed  that,  if  much  of  the  volume  of  the  chamber  were 
to  be  affected  by  ionisation  nue  to  the  discharge  rather  than  by  the  point 
source  presented  by  c  spark  plug,  the  low  flame  speed  would  not  have  to 
bo  critical. 

The  investigation  was  conducted  by  means  of  a  laboratory  pulse 
generator  which  discharges  its  high  energy  by  means  of  electrodes  mounted 
in  a  "bomb"  chamber  containing  controlled  mixtures  of  aertonia  and  air  at 
various  pressures  and  temperatures . 
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APPARATUS 

• 

A  "bomb"  chamber  with  internal  configuration  appruxim.it  ing  the 
combustion  chamber  at  TOC  of  thn  LOS-465  engine  was  dcsi  mod  and  built 
by  nine  tho  upper  portion  of  a  piston  from  that  engine  and  a  flat  plate 
to  simulate  the  cylinder  head  (see  Fig.  5).  Several  variants  of  this 
head  were  built,  to  allow  for  testing  of  different  electrode  configura¬ 
tions.  The  mixtures  were  introduced  Into  the  "bomb"  from  high-pressure 
tanks  by  means  of  pressure-regulation  devices  and  valve  arrangements, 
as  shown  In  Figure  I. 

The  initial  plan  called  for  In-house  design  and  construction  of 
the  high-energy  pulse  generator.  It  was  found,  however,  that  the 
Tobe-b'sutschmann  Company  had  }ust  begun  to  manufacture  a  unit  represent¬ 
ing  the  forefront  of  the  state  of  the  art.  Their  unit  is  capable  of 
producing  up  to  10  joules  at  up  to  200-fcV  output  pulse,  with  a  duration 
of  lO”^  seconds.  The  unit  renuires  only  a  10-kV,  low  current  supply. 

In  physical  size  (exclusive  of  power  supply/  it  measures  approximately 
8  inches  in  diamrter  by  4  Inches  in  height. 

A  su.udard  coil  and  breakn  ignition  system  was  used  for  com¬ 
parison.  A  Delco-Remy  128  -o>l  we*  ti*ed.  By  increasing  the  primary 
voltage,  output  energy  up  to  about  60  millijoules  was  made  available. 
Spark-plug  gaps  ranged  from  0.020  vo  0.060  Int'cs. 

Dynamic  pressure  measurements  were  taken  with  a  Kistler  natural- 
ouartz  piezoelectric  transducer,  an  amplifier  calibrator,  and  a 
Tektronix  oscilloscope.  Pressure-time  recordings  were  made  by  means 
of  an  oscilloscope  camera  (sea  Fig.  4  for  typical  example). 
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TEST  PROGRAM 

In  order  to  test  the  purser,  Atmospheric  discharges  with  varying 
arc  length  were  tried.  Spectacular  arcs  up  to  10  Inches  in  length  were 
drawn  (see  Fig.  2).  When  of  any  length  in  excess  of  about  4  inches,  the 
arc  is  not  confined  to  a  single  main  bolt,  but  rather  is  spread  to  many 
fine  tendrils  following  the  lines  of  the  electric  field. 

Automotive  spark  plugs,  suitably  modified,  were  selected  for  use  as 
electrodes.  Modification  entailed  the  cutting  away  of  the  ground 
electrode  and  a  portion  of  the  thread-,  in  order  to  provide  spacing 
adcouate  for  prevention  of  flash-over. 

A  series  of  "open  air"  observations  was  made,  with  the  electrode 
mounted  In  the  head  plates  to  determine  the  path  of  the  arc.  A  strong 
tendency  to  "short-cut"  through  the  head  plate  was  noted.  A  considerable 
measure  of  care  had  to  be  exercised  in  locating  the  electrodes,  in  order 
to  centralize  the  discharge  in  the  combustion  chamber  ard  provide  a  long 
discharge  path  without  causing  the  arc  to  short-cut  either  to  the  head 
or  piston  crown  (see  Fig.  3). 

After  a  f ew  firings  in  the  pressurized  chamber,  the  performance  of 
the  system  deteriorated  rapidly.  It  was  found  that  the  ceramic  insula¬ 
tion  at  the  spark-plug  nose  had  been  punched  through  by  the  discharge, 
reducing  the  arc’s  length  and  severely  masking  it.  After  some  trial  and 
error,  it  was  found  that  certain  series  of  Autoiitc  plugs  have  signifi¬ 
cantly  thicker  insulation  at  the  nose.  These  performed  satisfactorily 
throughout  the  remainder  of  the  test  program. 

An  effort  was  made  to  use  hydrocarbon  fuel,  in  the  form  of  gasoitne, 
for  reference  and  cal ibretion.  Since  the  test  equipment  was  designed  for 
gaseous  fuel,  difficulty  in  controlling  the  air-fuel  ratio  and  the  degree 
ol  vaporization  rendered  this  approach  Impractical .  The  obvious  re¬ 
course  was  to  use  a  gaseous  fuel,  thereby  permitting  familiarization  with 
thr  gas-handling  eauipment  for  air-fuel  ratio  control  («ee  Fig.  i). 

Thu*,  n  though  hydrocarbon  fuel  wr>uld  have  provided  better  reference 
data,  hydrogen  had  to  be  used  for  check-out  and  calibration. 
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The  optimum  air-fuel  ratio  with  hydrogen  was  determined  experimen¬ 
tally.  A  determination  was  made  of  the  effects  on  total  combustion  tine  of 
total  chamber  pressure  and  of  total  discharge  energy  at  lean  air-fuel 
rat  ios . 

Ammonia  combustion  was  then  attempted.  Again  the  optimum  air-fuel 
ratio  was  experimentally  determined.  The  dlfferaejce  in  total  combustion 
time  between  hydrogen  and  o<monla  fuel  was  very  clearly  evidenced,  not 
only  by  the  oscillograph  record,  but  by  observations  through  a  quartz 
window  in  the  combustion  chamber.  Reliable  awwnia  combustion  was 
achieved  with  the  high-energy  system,  for  all  pressures  in  the  range  from 
200  to  600  psig. 

The  effec  .>?  discharge  energy  at  constant  total  pressure  was  In¬ 
vestigated  next,  for  air-fuel  ratios  at  or  near  stoichiometric.  It  was 
found  that  combustion  time  did  not  increase  signif icent ly  even  at  the 
lowest  energy  of  which  the  pulse  generator  was  capable  (about  5  ioules). 

For  the  sake  of  convenience,  these  first  tests  were  performed  with 
the  bomb  and  gases  at  room  temperature.  It  was  felt  that  these  conditions 
(and  the  lack  of  turbulence  and  sw;H),  although  not  realistic,  were  very 
much  more  severe  than  those  in  an  actual  engine.  Chamber  temperature 
was  then  elevated  by  insertion  into  an  oven  —  with  a  practical  limit  of 
710  R  set  to  prevent  damage  to  equipment.  The  temperature  dependence  of 
combustion  time  and  pressure  rise  was  thus  determined. 

Up  to  this  time,  all  tests  had  been  performed  with  a  head  plate 
which  provided  a  3/16-inch  gap  3/8-lnch  below  the  head  (see  Figs.  5,  6, 
and  7).  Row  a  second  heed  plate  was  installed  In  which  the  arc  was  placed 
deeper  it.  t:_.?  chamber  (y  inch)  and  lengthened  to  3/8  Inch;  and  the  tests 
ware  repeated. 
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RESULTS 


Experimental  results  arc  plotted  in  Figures  8  through  15-  Figures 
8,  ?,  and  10  show  the  results  of  calibration  runs  with  hydrogen. 

Figure  II  indicates  the  determination  of  best  fuel-air  ratio,  on 
the  basis  of  minimized  combustion  time.  This  was  found  to  be  slightly 
in  excess  of  iS^percent  ammonia  by  pressure. 

The  effect  of  total  initial  pressure  on  the  total  combustion  time 
is  relatively  small  in  the  region  tested  (sac  Fig.  I?).  The  variation 
is  on  the  order  of  ±10  ml  1 1 (seconds  about  the  average  of  about  100  msec, 
showing  no  specific  trend. 

The  effect  of  discharge  energy  on  total  combustion  time,  shown  in 
Figure  13,  is  also  small  in  the  range  provided  by  the  pulse  generator. 

The  energy  provided  by  a  coil  and  breaker  ignition  system,  however,  was 
insufficient  to  initiate  combustion  at  all.  This  seems  to  indicate  that 
the  pulse  generator  is  much  more  productive  than  necessary  to  obtain 
reliable,  rapid  combustion.  Perhaps  the  most  Interesting  feature  of 
Figure  13,  however,  is  the  large  difference  in  total  combustion  time,  at 
all-  energy  values,  between  the  small-  and  large-gap  electrode  configura¬ 
tions.  The  combustion  time  with  small  gap  is  nearly  halved  in  going  to 
the  larger  gap.  This  sucoorts  the  initial'  premise  that  increasing  the 
amount  of  gas  affected  by  the  discharge  (by  increasing  arc  length)  would 
tend  to  counteract  the  slow  flame  speed  of  ammonia. 

The  effect  of  temperature  on  combustion  lime  is  plotted  in  Figure 
14.  Combustion  time  decreases  significantly  with  temperature,  for  both 
large  and  small  discharge  paths.  Of  more,  interest  than  the  actual  values 
In  the  range  tested  is  the  trend  indicated  by  the  curves.  Polytropic 
compression  (k»1.35)  with  compression  ratios  of  10  and  20  yields  com¬ 
pression  end  temperatures  of  about  1200°  and  1460^R  respectively.  A 
conservative  extrapolation  of  the  curves  in  Flnure  >4  yields  results  on 
the  order  of  a  few  milliseconds  —  l.e.,  the  practical  ranee  for  an  engine. 


Perhaps  the  most  interesting  combustion  characteristic  is  the  rate 
of  pressure  rise.  !n  Figure-  1 5 .  pressure  rise  per  rj  *  I > second  is  plotted 

against  temperature,  with  the  Initial  pressure  and  discharge  length  as 
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the  parameters.  The  pressure  rate  is  wry  strongly  affected  by  the  dis¬ 
charge  length,  the  longer  path  yielding  between  two  and  three  times  the 
rate  of  the  shorter  path.  Also,  the  increase  in  rate  with  temperature 
Is  quite  strong.  These  results  can  be  interpreted  In  terms  of  an  engine 
by  relating  pressure  rise  per  millisecond  to  crank-angle  degrees  of  an 
anginc.  At  2,000  rpm,  one  millisecond  Is  equivalent  to  1.2  degrees  of 
rotation.  In  a  diesel  engine,  pressure  rates  of  3  to  6  atm  (44  to  8€  psi) 
per  degree  are  normal.  An  engine  with  10:1  compression  ratio  will  have 
compression  end  conditions  of  1200°R  and  330  psig.  Likewise,  the  20:1 
engine  will  have  1460  R  and  840  psig.  Extrapolating  the  experimental 
data  conservatively,  and  divining  by  1.2  to  convert  to  eran4  angle  at 
2,000  rpm,  we  obtain  6  and  iO  atm  (88  and  147  psi)  per  degree  crankshaft, 
respectively,  for  the  two  compression  ratios. 

In  summary,  the  high-energy  system  produced  reliable  combustion 
under  much  less  favorable  conditions  than  one  would  expect  in  an  engine 
even  under  adverse  conditions  such  as  cold  starting,  and  —  in  fact  — 
even  under  circumstances  in  which  a  standard  ignition  system  failed  to 
produce  ignition  at  all.  Extrapolat ion  of  results  of  tests  indicates  that 
a  high-energy  system  can  reduce  total  combustion  time  to  appropriate 
values.  If  this  energy  is  rcinnsed  through  a  large  portion  of  the  mixture, 
pressure  rates  can  be  raised  sufficiently  for  engine  operation,  even  in 
a  test  chamber  without  turbulence,  and  at  low  temperatures. 
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CONCLUSIONS 

It  has  been  demonstrated  to  our  satisfaction  that  the  initial 
assumptions  (witn  rcgaiu  to  the  effect  uf  a  it«gn-«n*f gy  line  discharge) 
hold  •«  and  tht-t  aaaonia-conb-jtion  ch-ir  jclc'' •  , ore  er-hjoced  by  the 
use  of  a  high-energy  line  discharge. 

It  is  highly  significant  that  the  high-energy  system  ignited  the 
mixture  at  conditions  under  which  the  coll  and  breaker  system  failed 
completely.  It  is  also  very  significant  that  increasing  the  arc  length 
reduced  combustion  time  drastically. 

To  arguments  against  the  extrapolation  to  engine  temperatures, 
from  test  temperatures,  it  may  be  pointed  out  that  these  extrapolat Ions 
were  made  very  conservatively.  (Furthermore,  even  at  the  temperature  of 
7I0°R  the  pressure  rate  is  on  the  order  of  3  atm  per  millisecond,  and  the 
results  show  a  safety  margin.) 

There  are  strong  indications  that  the  energy  needed  to  yield  the 
improvement  in  combustion  is  less  than  that  produced  by  the  laboratory 
pulse  generator,  and  that  a  unit  for  engine  use  could  be  smaller.  In 
preliminary  discussions,  representatives  of  the  Tobe-Deutschmann  Company 
have  indicated  that  a  unit  of  reduced  energy  capacity  can  be  but  It  with 
repetition  rates  commensurate  with  engine  requi cements,  without  encounter 
with  any  great  technical  difficulty.  Furthermore,  other  manufacturers  of 
high-energy  pulse  systems  are  developing  units  which  may  be  applicable, 
anc-  available,  soon. 
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RECOMMENDATION 

in  the  light  of  the  promising  results  of  these  tests,  it  is  recom- 
mended  that  a  pulse  generator  with  attendant  eouipment,  suitable  for 
engine  operation,  be  assembled  and  tested.  Several  engines  should  be 
selected  for  the  purpose  and  modified  to  incorporate  electrode  location 
and  spacing  which  will  enhance,  to  the  maximum  possible  extent,  the 
properties  found  des i rable  In  this  study. 
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RADIO-FREQUENCY  DISSOCIATION  STUDY 

During  early  discussions  of  the  problems  presented  by  ammonia- 
fueled  internal  combustion  engines,  attention  was  centered  on  the 
question  of  ignition.  The  larger  ignition  energy  and  the  slower  flame 
speed  for  ammonia-air  as  compared  with  hydrocarbon-air  mixtures  suggested 
consideration  of  alternatives  to  the  conventional  spark  ignition.  One 
scheme  that  appeared  to  be  attractive  was  that  of  dissociating  a  part  of 
the  incoming  nitrogen  in  a  manifold  by  means  of  an  electrodeless  RF 
discharge.  The  atomic  nitrogen  would  then  be  distributed  to  the  engine 
cylinders  along  with  the  ammon'a-air  mixture.  Upon  recombining,  the 
nitrogen  dissociation  energy  would  be  recovered  and  would  be  available 
for  ignition  of  the  ammonia.  In  this  way  large  amounts  of  ignition  energy 
could  be  added  to  the  air-fuel  mixture  and,  most  importantly,  the  energy 
would  be  distributed  over  the  entire  volume  .of  reactants,  instead  of 
being  concentrated  in  a  small  region  of  space  around  the  spark  discharge. 

Experience  with  RF  dissociation  of  pure  nitrogen  indicated  that 
the  RF  energy  source  was  not  particularly  difficult  to  construct,  that 
the  atomic  nitrogen  has  a  sufficiently  long  lifetime,  and  that  the 
recombination  can  be  easily  stimulated  by  interposing  a  refractory  surface 
into  the  dissociated  nitrenen  stream.  However,  it  is  important  to  note 
that  in  the  RF  discharge- region  relatively  high  temperatures  exist;  and 
this  has  the  effect  of  prolonging  the  atomic  nitrogen  lifetime.  In 
view  of  the  differences  between  the  RF  dissociation  of  pure  nitrogen 
and  the  use  contemplated  here  as  an  aid  to  the  ignition  of  an  ammonia- 
air  mixture,  it  was  felt  that  theoretical  calculations  should  be  undertaken 
before  even  consideration  of  any  practical  implementation  of  the  idea. 

Within  the  limitations  imposed,  a  modest  program  was  undertaken 
focusing  on  two  basic  problems: 

a.  The  time  scale  of  the  nitrogen  recombination  as  a  function  of 
temperature  and  initial  atomic  nitrogen  density 

b.  Significant  effect  of  the  presence  of  oxygen  on  the  time  scaie 
of  the  recomb inat ion 
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These  questions  have  been  examined  by  solving  appropriate  sets  of  differ¬ 
ential  equations  that  express  the  rates  at  which  the  various  chemical 
reactions  proceed.  By  carrying  out  such  calculations  parametrically,  in 
terms  of  the  degree  of  dissociation  and  temperature,  considerable  insight 
can  be  obtained  into  the  processes  taking  place.  In  Section  A  the 
kinetics  of  the  dissociated-nitrogen  air  system  are  discussed.  This 
includes  all  of  the  possibly  important  reactions  and  their  rate  coeffi¬ 
cients.  In  Section  B  a  number  of  numerical  results  are  presented  for 
a  series  of  increasingly  sophisticated  models  used  to  describe  the  time 
behavior  of  the  dissociated-nitrogen  air  system.  Conclusions  and  recom¬ 
mendations  are  contained  in  Section  C.  Finally,  in  the  Appendix  of  this 
report,  some  of  the  numerical  techniques  employed  are  briefly  summarized. 
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KINETICS  OF  THE  01SSQCIATED-NITR0GEN  AIR  SYSTEM 

It  is  assumed  that,  by  means  of  an  RF  discharge  operating  in  air, 

N„  is  dissociated  to  atomic  nitrogen. 

Na  -*  N  +  N  (I) 

Furthermore,  it  is  assumed  that  the  corresponding  oxygen  dissociation 
does  not  proceed  significantly  until  al>  of  the  molecular  nitrogen  is 
dissociated.  In  addition,  we  neglect  the  following: 

a.  Mass  diffusion  from  the  manifold-reactor  region 

b.  Wall  effects 

c.  Heat  addition  due  to  chemical  reactions 

Thus,  we  have  as  the  initial  state  of  the  manifold  a  mixture  of  N,  ,  0;,  • 
and  N  .  The  pressure  is  likely  to  be  atmospheric  and  the  temperature 
*  30CPK  (5*»0°R)  depending  on  the  degree  of  heating  resulting  from  the  R» 
discharge. 

Primary  Reactions 

The  lifetime  of  the  atomic  nitrogen  is  of  major  interest.  We  can 
approach  this  problem  most  simply  if  we  restrict  ourselves  to  considera¬ 
tion  of  the  reactions  that  can  be  supported  by  the  initial  mixture  only. 

One  class  of  such  reactions  is  nitrogen  recombination  whit!-'  usually 
proceeds  via  the  three-body  process 

N  4  N  +  M  -  Na  +  M  (2) 

where  M  can  be  H  ,  H5  ,  or  05  .  A  very  complete  review  of  the 
reaction-rate  literature  relating  to  this  process  can  be  found  in  Cary. 
Reaction  rates  when  the  third  body  is  N  or  N..  are  well  known;  the  reac¬ 
tion  rate  when  the  third  body  is  0.,  has  not  been  specifically  investi¬ 
gated,  but  one  can  assume  that  it  is  similar  to  the  rate  when  the  third 
body  is  Na  .  The  rate  of  the  radiative  recombination  process 
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N  +  N  -  N;!  +  hv  (3) 

is  negl  lgl!jle.^ 

Atomic  nitrogen  can  be  lost  net  only  by  nitrogen  recombination  but 
also  by  oxidation.  This  proceeds  either  to  NO  with  atomic  oxyftn  as  a 
product , 

N  +  0g  -*•  NO  +  0  W 

or  radiativ«')y  to  NO.,  , 

N  +  0g  -*  N0?  +  hv  15) 

The  rate  for  (4)  has  been  investigated  by  a  number  of  workers  and  can  be 
considered  to  be  reasonably  well  known. ^5*28,32  rata  tie 

4 

radiative  process  (5)  has  been  estimated  by  8ortner. 

Secondary  Reactions 

Note  that  (4)  and  (5)  form  0,  NO  ,  and  NO.,.  If  the  degra*  of 
nitrogen  dissociation  is  small  and/or  the  net  effect  of  (2)-(5)  does  not 
result  in  a  significant  ioss  of  N  then  the  formation  of  0  ,  NO  ,  and 
NQg  can  be  ignored.  However,  if  significant  amounts  of  C  ,  NO  ,  and 
N0a  are  formed,  then  their  reactions  must  be  considered  also.  This 
involves  a  large  number  of  reactions,  potentially,  although  one  expects 
to  find  only  a  few  that  are  actually  dominant.  The  secondary  reactions 
can  be  divided  into  these  classes:  ; 

a.  Reactions  of  0.  or  Ng  (which  are  expected  to  be  present 
in  substantia)  amounts)  with  0  ,  NO  ,  or  N0? 

b.  Reactions  requiring  N  (which  may  not  be  present  in  a  sub¬ 
stantial  amount)  and  0  ,  NO  ,  or  N02 

c.  Reactions  reouiring  0  ,  NO  ,  and  N0^  only 

d.  Reactions  requiring  0,  and  Na0  ,  which  are  formed  as 
products  of  reactions  in  classes  a  and  b 
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These  are  now  considered  in  detail. 


a.  Reactions  of  02  or  N_  with  0  ,  NO  ,  or  NO 


Tbe  three-body  reaction 


NO  +  NO  ♦  0-  -  NO-  +  N0„ 


ha*. been  investigated  by  Harteck  and  Dondes.  °  Ozone  formation  proceeds 
ei tiler  by  the  three-body  process 


or  relatively, 


o  +  o2+m-o3  +  m 


f  +  0-  -*  Oj.  +  hv 


1  39 

The  former  has  been  studied  by  Axworthy  and  Benson,  Zaslowsky  et  al, 

and  Campbell  and  Nudelman.^  The  reaction  rate  for  the  radiative  process 

L 

has  been  estimated  by  Bortner.  Finally,  there  are  the  oxidation  processes. 


0  +  N  -♦  NO  +  N 


0  +  Nc  -*  Ns0  +  hv 
0  +  NB  +  M  -  Ng0  +  M 


The  two-body  rearrangement  reaction  (9)  has  been  extensively  studied  by 

2 L  26  37 

Kaufman  and  Kelso,  Kistiakowsky  and  Volpi,  Wray  and  Teare,  Clyne  and 

Th'vsh,^  Phillips  and  Schiff.^  The  two-body  radiative  rata  for  (10)  has 

been  estimated  by  Bortner.1*  The  three-body  process  (It)  has  been  studied 

by  Hartecic  and  Dondes.1^ 


-;t*y 

b.  Reactions  of  N  with  0  ,  N0% 

Tha  leaction  of  N  and  0  to  M  ten  take  place  either  by  a  three- 
body  process, 
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(12) 


or  radiatively. 


N  +  0  -  NO  +  hv  (13) 

The  rate  of  (12)  has  been  studied  by  Harteck,  Reeves,  and  Hannella,^ 

2  21 

Barth,  and  Kaplan  and  8arth  while  the  rate  of  (13)  has  been  estimated 
by  Bortner. 

The  reaction  of  N  and  NO  can  take  place  either  by  the  rearrange¬ 
ment  process 


N  +  NO  -  N.  +0 

(14) 

studied  by  Kistiakowsky  and  Volpi,2^  Kicolet  and  Alkin 
Herron,'^  Clyne  and  Thrush,^  and  Phillips  and  Schiff^ 

32 

,  Kaufman  and  Kelso 

-  or  radiat ively. 

N  +  NO  -  rik0  +  hv 

05) 

for  which  the  rate  has  been  estimated  by  Bortner. ^ 

The  reaction  of  N  and  N0Z  can  proceed  to  a  variety  of  products: 

N  +  NO..  -  NO  +  NO 

(16) 

N  +  NO;,  -  NgO  +  0 

(17) 

N  +  NIL  -  Ns  +  0  +  0 

(18) 

N  +  NO,  -  N.  +  0, 

C  3 

(19) 

Reactions  (16)-(18)  have  been  studied  by  Harteck  and  Dondes^*'®  while 
the  rate  of  (19)  has  been  estimated  by  Bortner.^ 


c.  Reactions  reauiring  0  ,  NO  ,  NQ?  only 

Oxygen  recombination  can  take  place  either  by  a  three-body  process  , 

0  +  0  *  M  0-  +  M 


(20) 
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or  radiatively, 

0  +  0  -  0;  +  hv  (21) 

14 

The  three-body  process  has  been  studied  by  Golden  and  Myerson,  Kaufman 

47  34  29 

and  Kelso,  Reeves,  Mannella,  and  Harteck,  Morgan,  Elias,  and  Schiff, 

'  21 

and  Kaplan  and  Barth.  The  radiative  process  has  been  considered  by 
Nicolet^  and  Ni colei  and  Mange. 

0  and  NO  react  In  these  three  different  ways: 

0  +  NO  +  M  -  N02  +  M  (22) 

0  +  NO  -  N0S  +  hv  (23) 

0  +  NO  -  0.  +  N  (24) 

24 

The  three-body  reaction  (22)  has  been  studied  by  Kaufman  and  Kelso, 

2-1  1 2 

Kaufman,  and  Ford,  Doyle  and  Endow.  The  radiative  path  to  N0S  ,  (23), 

has  been  treated  by  Ciyne  and  Thrush, 10  while  the  rearrangement  reaction 

24  37 

(24)  has  been  examined  by  Kaufman  and  Kelso,  Wray  and  Teare,  and 
Kist’akowsky  and  Volpi.^"* 

Finally,  0  and  N0?  can  react  by 

N0?  +  0  ~  NO  +  0..  (25) 

12 

Reaction  rates  have  been  reported  by  Ford,  Doyle,  and  Endow,  Harteck  and 
Dondes,^  Kaufman  and  Kelso, ^  and  Phillips  and  Schiff. ^ 

d.  Reactions  requiring  0.  and  N  0 

The  two  reactions  of  interest  that  consume  03  are 

0  ♦  0,  -  0.  +  0.. 

«  »  , 

NO  4  C3  -  NO.  •»  Ca 

3 

The  first  of  these,  (26),  has  been  studied  by  Benson  and  Axworthy, 

Zar lowsky  et  al,^  Phillips  and  Schiff,^  and  Campbell  and  Nudelman,^ 


(26) 

(27) 
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The  second  reaction,  (27),  has  been  examined  by  Johnston  and  Crosby, 
and  by  Greaves  and  Garvin.^ 

Finally,  Nc0  can  react  with  0  in  the  following  two  possible  ways: 

0  +  Nj.0  0.  +  N?  (28) 

0  +  Na0  -  NO  +  NO  (29) 

The  first  of  these,  (28),  has  been  studied  by  Bradley  and  Kist iakowsky ,** 

1 1  22 

Fenimore,  and  Kaufman,  Gerri,  and  Bowman.  The  second  reaction,  (29), 

22 

has  been  considered  by  Wray  and  Teare  '  and  by  Kaufman,  Gerri,  and  Bowman. 


Other  Reactions 

There  are  a  number  of  other  possible  reactions  that  are  not  believed 
to  be  important  at  normal  temperatures  and  pressures.  While  it  is  con¬ 
ceivable  that  they  may  occur  during  the  high-temperature  and  high-pressure 
conditions  attendant  upon  the  ammonia  combustion,  the  low  concentrations 
of  the  reactants  probacy  still  justifies  their  neglect.  These  reactions 
are 


N0+H-N  +  0  +  M 

(30) 

N03  +  H  -  0  +  NO  +  M 

(30 

0,  +  M  0  +  0.  +  M 

(32) 

N.,0  +  M-  0  +  N  +  M 

(33) 

37  13 

Wray  and  Teare  '  and  Freedman  and  Oaiber  have  considered  (30)  while 

C  1  J  || 

Bradley  and  Kistiakowsky'7  and  Fenimore  have  considered  (3 1).  Bortner 
has  commented  on  (32)  and  (33) . 

Finally,  while  coilisional  dissociation  of  N2  and  0a  , 


N3  +  K-  N  +  N  +  M 
0a  +  M-  0  +  0  +  M 


(34) 

(35) 


obviously  does  not  proceed  rapidly  at  normal  temperature  and  pressure, 

8 

these  react  ions  may  play  a  role  during  combustion.  The  work  of  Cary 
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in  connection  with  (3*0  has  already  been  mentioned.  Studies  of  the 

35 

collisions!  dissociation  of  0  have  been  ma^e  by  Rink,  Knight,  and  Duff, 

6  2  y  *  3A  33 

Byron,  Matthews,  Schexnayder  and  Evans,  and  Wray. 

Reaction  Rates 

Before  quantitative  calculations  can  be  made  of  the  nitrogen  recom¬ 
bination  times  resulting  from  the  preceding  reactions,  one  must  have 
numerical  values  for  the  reaction  rates.  These  have  been  obtained  from 
several  secondary  sources  such  as  published  review  papers  and  unpublished 
industrial  reports.  These  in  turn  reference  the  primary  publication  of 
the  data.  The  reaction  rates  are  presented  below  in  tabular  form,  along 
with  references  to  both  the  primary  .nd  the  secondary  source  of  the 
information.  In  some  cases  the  information  is  available  as  a  function 
of  temperature,  while  in  others  a  single  value  for  a  specific  temperature 
is  given.  In  some  cases  rather  considerable  differences  in  the  reaction 
rates  may  be  observed.  These  differences  serve  to  emphasize  the  uncertainty 
attached  to  this  type  of  information,  sometimes  of  the  order  of  several 
orders  of  magnitude.  The  physical  origin  of  these  differences  lies  in 
the  difficulty  of  ascertaining  in  a  given  experiment  the  degree  of 
electronic  excitation  of  the  reactants  end  their  velocity  distribution. 

Thus  different  experimental  measu''emerts  of  what  is  presumed  to  be  the 
same  physical  quantity  sometimes  yield  rather  dissimilar  results.  In 
such  a  case,  however,  it  is  often  possible  —  by  judicious  choice  of  values 
for  a  calculation  —  to  set  upper  or  lower  limits  for  a  particular  reaction 
rote  and  thereby  establish  its  importance  cr  unimportance. 

Reaction-rate  tables  follow. 


18 
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TABLE  I 

REACTION  RATES  FOR  THE  OISSOC IATED-NI  TROr.EN  AIR  SYSTEM 
(Primary  R«*acl  ions) 


P*mrt  Inn 

tqiMt  Inn 

Nlimlin 

T  mp . 

Pr  !m.v  y 
Rrfplpiirr 

Snrnnrfn'y 

#«>frirnrn 

N4N4M-N..4N 

(M-N.) 

2 

»,«IO*,n'n-'i  T*'  tm'/src. 

- 

»is.i  ;.'»2 

4 

(M-N) 

S«»0*J°,n  T*» 

- 

- 

4 

(H/N.N,) 

j.tO-JO.05  T-l 

- 

- 

it 

1.4*10*” 

300"  it 
(S'-OA) 

7 

2k  10*” 

300^  It 

iW*) 

45 

'•3 

N«N  N.*hv 

3 

- 

mm 

- 

4 

K«0.*0*0 

it 

33.10*” 

- 

25 

2 

1  5*10-”  r*6?00/"T 

— 

2S.3? 

«*S 

Ti/»  _ 


24.2* 


4 
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B 


TABLE  7 

REACTION  RATES  FOR  THE  DISSOC IATED-NITROGEN  AIR  SYSTEM 
(Secondary  Reactions  -  Type-  a) 


React Ion 

EqiMt  Inn 
Number 

Rate 

Temp. 

P»  Invtry 
Reference 

Secondary 

Reference 

N0tN0+0.  -NO^NO, 

6 

lft-37'0.7  6. 

10  i.m  /sre 

300  K 
(540"*) 

18 

2*3 

0+0.,+M  -0,+M 

7 

K**3'* e.‘/w 

300'‘K 

(S'*  on*) 

1 

«*3 

6.8..0*31*  e300/T 

- 

3 

2 

3.10*  3,1 

- 

3.29.7 

A 

(HO^  -o4*h. 

e 

I0*20‘J  c<n3/see 

- 

- 

U 

0*N;>  -N0*N 

s 

,0-IO*C.5  r-37S00/T 

cm3 /set 

- 

2**. 26.37.9.33 

(* 

0*N,  -NyO+hv 

10 

I0-J0'*  r.J/sec 

- 

- 

i* 

0*Hrf+K  -H  O+M 

II 

7X.O’33  e*,200/RT  t 

_ 

- 

16 

i* 

rm  /sec 
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TABLE  2  (ContM) 

REACTION  RATES  FOR  THE  DISSOClATED-NITROGEN  AIR  SYSTEM 
(Secondary  Reactions  -  Type  b) 


React  Ion 

Equal  Ion 
Number 

Rale 

Temp. 

Pr tmary 
Reference 

Secondary 

Reference 

N40*H  •NO+H 

12 

SxlO'3"0^  T',/2  cm6/*ee 

- 

17 

*.2 

!.5*tO*32 

300°  K 
(5W'R) 

21 

2 

N+0  4f0*h>- 

13 

10*2J12  cm3/sec 

- 

- 

* 

tMNO 

1* 

2.2xl0*,,‘°’5  cm3 /sec 

- 

2*. 19.9.33 

* 

S.JxlO*"  , 

JOO^K 

i5-«e’R) 

26 

2 

1.5*10*"  T»'? 

- 

76.32 

*3 

H*H0 

15 

IO*20'2  cm3/sec 

- 

- 

U 

N4N0_,  -NO»NO 

16 

5..0*''“0-5  cm3/sec 

300°  K 
(5*0  *) 

18 

*3 

N4N0,  •N.,0»0 

17 

J.2xl0*',l,°  ^  cro3/sec 

300°  K 
(5*0'R) 

18 

*3 

N4M0_  -H.,4040 

18 

2x10*,MO'5  cm3/sec 

300°  K 
(5*0"R) 

18 

*3 

N4N0>,-N..40j> 

19 

2x10****"  cm3/sec 

- 

- 

A 

21 
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TABLE  2  (Cont»d) 

REACTION  RATES  FOR  THE  OISSOCIATEO-NITROGEN  AIR  SYSTEM 
{Secondary  Reactions  -  Type  c) 


Reaction 

tqn.it  Ion 
Numt.er 

•  .111’ 

Temp. 

f  r  inwry 
Reference 

Secondary 

Reference 

0+0+M  •<>,♦*  (M-0.) 

20 

6x,0-37*0.5  t-I/2  ef/w 

- 

JM7.3fc.79 

4 

(M-0) 

2xIO-3I»0.5t-»/? 

- 

- 

4 

(M/0.0,) 

|0-32*0.5  t-I/2 

- 

- 

4 

1.6*»0'37 

0 

jorTK 

CcWr) 

29.21 

? 

3*I0*33 

300°  K 
(SfcO'R) 

14 

43 

2»0  -0^4hv 

71 

»o*27'7 

- 

31,30.48 

4.43 

(M  M... A)  22 

5. 5*  10* 37  cm* /sec 

300"k 

(r.fc0'A) 

23 

43 

5*I0'37 

- 

24 

4 

5.?«IO*32 

300°K 

I5'*0'A) 

12 

2 

0*H0  •*0.,4hv 

2J 

IO*20*7  cm3 /sec 

- 

10 

4 

04N0  ‘O^+H 

24 

)0-l3‘0.5  t»/2  e-3B600/AT 

«3/i« 

- 

24.37.25 

4 

NO<,40<NO»0.> 

25 

10*  *2  c">3/»ec 

300"  K 
(5'«0"R) 

18 

43 

r-h?2 
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REACTION  RATES  FOR  THE  DISSOCIATED-NITROGEN  AIR  SYSTEM 


(Secondary  React 

ions  -  Type 

d) 

React  Ion 

tqti.it  Ion 
Numher 

Rate 

Temp. 

Primary 

Reference 

Secondary 

Reference 

040,  0,40, 

26 

5W"  e-3000/T  rmJAee 

- 

3 

2 

7.I.IO',,>0-5  e-,#100/T 

- 

3.39.33.7 

4 

N0+0, '410,40, 

27 

8„.0*,J  T,/2  e*2500/W, 

«r/$fc 

l98-7)0°lt 
(356-4  |4°R) 

20.15 

43 

04M.P-0,4H, 

28 

^.O-.MI  ,-27000/RT 

C»7»fE 

- 

5.22.11 

4 

04R,0  “04N0 

29 

'  ..7KlO*,0“  e*28000/" 

- 

37.22 

4 

e*J/s« 
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TABLE  3 


UNIMPORTANT 

REACTIONS  FOR  THE 

DISSOCIATEO-N ITROGEN 

AIR  SYSTEM 

Reaction 

Equat ion 
Number 

Primary 

Reference 

Secondary 

Reference 

N0+M-N+0+M 

30 

37,13 

4 

NO.j+M-0+NQ+M 

31 

5,  H 

4 

0  +M-0+0_+M 

32 

4 

N..0+M-0+Na+M 

33 

- 

N^+M-N+N+M 

34 

40,41,42,8 

4 

0  +M-0+0+M 

35 

35,6,27,36,38 

4 

i 

i 
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CALCULATION  OF  NITROGEN  RECOMBINATION  TIMES 

We  hav»-  considered  four  models  for  the  chemical  reactions  occurring 
after  a  given  amount  of  nitrogen  is  dissociated.  The  first  two  models 
lend  themselves  to  analytic  solutions;  the  last  two  were  solved  numerically. 

Model  1 

In  Model  1  we  consider  the  recombination  of  N  to  N.  by  a  three- 

d  ' 

body  reaction  where  the  third  body  M  may  be  0;,  ,  N.,  ,  or  N  .  Furthermore, 
we  assume  that  the  concentration  of  the  third  body  is  constant. 

The  chemical  equation  is 

N  +  N  +  M-*N„+M  (36) 

and  the  associated  rate  equation  is 

-  k  n[M]  'infN)/  (37) 

where  n[X]  Is  the  concentration  of  X  ,  (X  *  N  ,  N..  or  0.)  and  k  is 

c 

•the  rate  constant  for  this  reaction. 

if  we  dissociate  t  percent  of  the  Na  ,  then  the  initial  amount 

o:  N  is  given  by  2f/l00  times  the  initial  concentration  of  N,  .  The 

initial  concentration  of  N..  can  be  computed  by  noting  that  at  STP  there 
23  * 

are  6.02  x  10  molecules  of  air  per  22.4  liters.  Hence  the  concentrat ion 
of  air  is  given  by 

23 

A«P.1,,.X,  19—.  >  2.69  x  10*^  molecules  of  air  per  cm^ 

22.4  x  I0J 

If  we  assume  air  consists  of  20-percent  0-  and  80-perccnt  N  ,  the  initial 
concentration  of  N  ,  nQfN]  ,  is  given  by 

n  IN]  »  4.3f  x  10^  cm’^ 
o 


Here  we  have  assumed  standard  temperature,  T  ■  300o,r.  .  If  we  assume 
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the  concentration  is  Inversely  proportional  to  the  temperature,  then 

n  [N]  Is  given  by 
o 

nQ[N]  -  4.3  f  x  I017  (300/T) 


If  we  assume  nfM)  is  e  constant,  then 


n[M]  «  (2.69  x  I019)  $jr) 


The  rate  constant  k  is  given  by 


k  *  5.0  x  10~9®/T  cm^/sec 


Solving  equation  (37)  we  have 

n  In) 


nom 


1  +  n  fN]  n  k  t 
o 


(38) 


Inserting  numerical  values, 


n(N)  «  T1 


1.29  x  10 


20 


J-  +  (5-2  *  l°'2)  t 


T  ^  .j2 

Therefore,  for  large  values  of  time  (t  >:>  x  10  )  we  have 

b.Zr 


n IN]  =  T3  (2.48  x  107)/t 


Hence  the  amount  N  ,  ir  this  approximation,  is  essentially  Independent 
of  the  initial  amount  of  N  ,  decreases  as  l/t  ,  and  is  proportional  to 
the  square  of  the  temperature  (the  temperature  once  set,  Is  assumed 
constant) . 


Model  2 

In  Model  1  we  assumed  that  the  concentration  of  the  third  body 
remained  fixed;  in  fact  this  is  not  the  case,  for  each  time  you  form  an 


26 
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Ny  molecule  you  lose  ti«  atoms  for  each  molecule  gained.  Also,  there  Is, 
a  small  difference  between  the  rate  constant  when  the  third  body  Is 
or  N  ,  and  0^  . .  Model  2  takes  these  effects  Into  consideration. 

The  chemical  equation  is  the  same  as  before,  but  the  rate  equation 
is 

-  <kni>tNl  +  k„n  f»a)  ♦  k0n  (0al)  n[»)  (39) 


where  k  is  the  rate  constant  for  N_  or  ti  as  the  third  body  end  k 
n  p  o 

is  the  rate  constant  for  0p  as  the  third  body. 

The  solution  for  (39)  is 


where 


nlK)  -  ' 
B  -  not«l.) 


"°|M|  L* 

1+n  (N)B  K  t  I 
o  nu 

+  *  nofNl  +  k^  no{°*; 

k  *  5  x  l<f30/T 

n 

k  »  3  x  1 0"30/T 

o 


(40) 


& 


& 


& 


Therefore  as  the  concentration  decreases  we  have 


and  hence  the  solution  takes  the  same  form  as  (38)  with  nfM)  replaced 
by  B  . 

If  we  plotted  the  Mode)  I  or  Model  2  solution  on  a  logarithmic  scale, 
we  would  expect  asymptotically  to  have  a  straight  tine  with  slope  'I 
Figures  16  and  17  show  the  results  for  Models  «  and  2  for  10-percent  dis¬ 
sociation  at  T  *  300°K  and  at  T  «  600' K  .  Several  points  are  apparent 
from  these  figures.  First,  the  time  scale  for  the  recoT&inei ion  is  short. 
At  about  0.)  msec  the  atomic  nitrogen  concentration  has  decreased  by  two 


2/ 


V^Y^'l "  ’**  vfcjsg :  r "  _*  V-l*-  jg^Sw—s- 


. .  „  ;-*•'  3r-^?  :-  --  --^~  -  -==-  -  --*<• 

.-  -saJ-Crffr 


z£ 
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orders  of  magnitude.  By  this  t»m<*  the  asymptotic  behavior  has  beyn  achieved' 
so  that  nfNl  l/t  ;  hence  the  atonic  nitrogen  concentrat Ion  decreases  by 
another  order  of  magnitude  by  I  msec.  To  put  the  time  scale  into  proper 
perspective,  one  can  assume  an  average  linear  flow  velocity  from  the 
manifold  into  the  cylinder  of  about  100  m/sec.  If  the  distance  from  the 
atomic  nitrogen  source  to  the  cylinder  is  5  cm,  then  the  transit  time  is 
about  0.5  msec  and  so  a  sufficient  amount  of  atomic  nitrogen  must  sur¬ 
vive  to  times  of  this  order  of  magnitude  if  it  is  to  assist  in  the  ammonia 
combustion. 

The  second  major  point  is  that  high  temperatures  lead  to  higher  con¬ 
centrations  at  late  times.  This  is  due  to  the  t/T  dependence  of  the 
recombination  constant.  At  higher  temperatures  k  decreases  and  so 
the  recombination  proceeds  more  slowly.  The  manifold  temperature  will  be 
of  the  order  of  room  temperature  (TOOK)  while  the  cylinder  temperature 
at  peak  compression  may  be  of  700-80oPk.  It  was  with  these  temperatures 
in  mind  that  30(fK  and  a  conservative  600r K  was  assumed.  However,  the 
major  significance  of  the  600r'K  calculations  is  to  demonstrate  the  direc¬ 
tion  of  the  temperature  dependence  of  the  reaction  rates,  since  it  is  the 
lower  temperature  manifold  conditions  that  arc  expected  to  control  the 
atomic  nitrogen  recombination. 

The  third  point  is  that  aside  from  seme  minor  differences  for  times 

less  than  0.1  msec,  Models  I  and  2  qlve  identical  results  at  both  300° 

and  600°K.  This  similarity  is  whal  one  expects  when  the  degree  of 

dissociation  is  small.  If  n  fN]  is  not  larne,  the  third-body  concen- 

o 

tration  is  almost  constant.  Regardless  of  the  Initial  dissociation,  there 
will  always  be  some  late  time  when  nfNl  has  decreased  sufficiently  to 
a  point  where  the  third-body  concentration  is  no  longer  changing  signifi¬ 
cantly. 

Model  3 

Up  to  now  we  have  neglected  any  reactions  involving  0,  with  N  . 

The  two  most  significant  are 


N  +  0.  **  NO  ♦  0  k,  »  (3.3  x  l(f ,2)  e*3100/T  cm3 /sec 

N  ♦  0.  ~  HO.  +  hv  k.. 


I0*20  cm3/sec 


(41) 

(42) 


II 

0 


0 


Q 

0 


A 


* 

I 


D 


I) 

a 


28 


fi 

0 

fi 

fi 

0 

fi 

0 

c 
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Combining  these  reactions  with  (36)  We  are  led  to  the  following  set  of 
differential  equations: 


*  -  {knr»fNl  +  knnfiy+  koM'  f  nfNj;  ‘  (43) 

-  (k,  +  It,)  nfH)  nfO.  ] 

*  •»  »’ 

dnfO.) 

— ~~  -  -  (kj  +  k.  )  nfH)  nfOj,)  (44) 

dn  fN  )  •  , 

——  -  \kn  In)  +  k  nfH;,)+  k  M»  nfN]  (45) 

dt  *■  n  n  o  j  <- 


where 

M'  ■  nf0_l  +  nfNJ]  +  nfO)  +  n I N0„ ) 

?k  +k 

M'  -  nfCL.1  +  |?(nofN2J  -  "W)  +  n^N)  -  nlN] 


0 

B 

0 

8 

0 

0 

n 

V.  - 
0 

fi 


Equations  (43)  to  (45)  form  a  set  of  nonlinear  coupled  differential 
equations.  With  the  aid  of  a  subroutine  PRECQR  that  uses  a  predictor- 
corrector  technique  to  solve  differential  equations  (discussed  in  the 
appendix),  numerical  solutions  for  this  model  were  obtained. 

Figures  18  through  20  show  solutions  plotted  on  a  logarithmic  scale 
for  T  -  300°  K  and  60<f  K  with  initial  nitrogen  dissociation  of  I,  5.  and 
10  percent.  Note  that  there  is  an  appreciable  change  in  the  behavior  of 
the  solutions  for  different  values  of  temperature.  This  is  explained 
by  the  increase  In  the  rate  constant  for  reaction  (4 1 )  as  the  temperature 
increases.  At  T  *  300°K  ,  k:  *  1.07  x  10"^  ,  while  at  T  «  bOO^K  , 
kj  -  1.88  x  10* 14  . 

If  one  assumes  that  n(02)  is  roughly  constant  and  that  the  only 
reaction  of  any  consequence  is  (41),  then 


-k  p i 0  1  f 

=  rt  !  N  j  *> 

M 


n  (N) 


>  c- 


- .  -r  =  ---:=u._  ~  ~_-_<f^r-?*  ^1^31..,'“  r.-?^^v_^-^«.i*«  -  «e^*>5vL.  .. 
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Taking  the  togarithm  of  both  sides  we  have 


-  kjnlOjJ  t 


Therefore,  if  in  fact  equation  (4|)  does  dominate  the  reaction,  the  plot 
of  n[N]  vs.  t  on  semi-log  paper  should  be  a  straight  line  with  slope 
-(k,n{a.J)  ,  Figure  21  is  such  a  plot  for  1-percent  dissociation.  Its 

slope  is  -4.8  x  10  while  -^nfO.  1  it  approximately  -4.78  x  I0*\  Hence 
we  conclude  that  at  higher  temperatures  the  oxidation  reaction  (4|)  domin¬ 
ates.  At  the  lower  temperature  the  three-body  nitrogen  rccombi nation 
reaction  continues  to  dominate,  as  is  indicated  by  the  nearly  straight 
line  on  the  logarithmic  plot  that  is  identical  to  the  Model  2  result  in 
Figure  17  for  t  >  0. I  ms?c.  Reaction  (42)  has  a  rate  constant  that  is  too 
small  to  be  of  any  significance.  The  essentially  new  fact  revealed  by 
these  calculations  is  the  rapid  decrease  in  the  atomic  nitrogen  concentra¬ 
tion  for  t  -0.1  msec  at  higher  temperatures  due  to  the  oxidation  reaction 
(41). 


Model  4 

Model  3  took  into  account  the  major  reactions  concerning  N  ,  0a  , 
and  N.;  ;  in  Model  4  we  take  into  account  reactions  of  these  gases  with 
the  molecules  formed  by  the  reactions  of  Model  3.  The  chemical  equations 
are  the  three  from  Model  3  plus  six  reactions  involving  NO  ,  0  ,  NO,  . 

They  are 


N+N-HHtj+M 


N+C-  -NO+O 


N+0.  -*N0.,+h\ 


N.  +0+M-N.O+M 


N~  •KMNOtN 


kn  «  5.0  x  10~30/T 
kQ  *  3.0  x  10"30/T 
k;  .  3.3  x  10'12  e'3,00/T 


k.  «  10**u 

k  -  ?,0  »  iO-33  -*6060/T 


k.,  -  10‘,(J  e”37500/T 


30 


I 

I 

6 

0 

C 

e 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

n 


CAE  Report  Ho.  1054 
'  :n:*  Appendix  IV 

Volume  II 

R.-HMi.O+hv  ^  »  I0~20 

0Jr+0+H-,0rj+M  k,  .  3. 162  x  10*35 


0^  +0-0.  +hv 

0.  +N0+N0-N0  +N0 


k,  -  10 


ks-  10 


-20 


•37 


The  estociated  ••efe  eouet ion«  arm 


dnlN) 

dt 

dnfN.  ] 
dt 


“  Ao  *  +  A< 

Aq  “  ^  "  A4  A, 


dnlOgl 

dt 


A  “A..  -  A.  -A  -  A 

1  *•  '  7  g 


dnfMQ] 

dt 


A;  +  A4  -  Aa 


dnfNO.,  J 
dt 


A  +  2A. 


dnfO] 

dt 


A  -  A,  -  A  -  A  -  A  -  A 


dnltf.O] 

dt* 


A  +  AB 


^  -  A  ..A, 


where 

Ao  “  ikn^nfNl  +  n  iN„  1  >  +  ko(H-nfN]  -  nfN.J);  ( .  [  0 )" 
A,  »  k,n[N]  ntaj 
Aa  »  kj,n{N]  nlOy) 

Ag  **  k,*H  •n[Ns].n(0] 


I 
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Ai  *  k^nf^J.nCOj 
A,  »  k»  n  (N,  ]  nfO] 

A,  ■  kp.nlC., J  n(0)H 
k,  ■  k^nlOyJ  nlOj 
A8  -  k8n(0j  (n(H0]f 

end  M  ■  summation  of  the  concentrations  of  ell  the  constituents. 

This  set  of  differential  equations  was  solved  numerically  using 
PRECOR.  Figure  22  compares  Model  3  and  Model  4  at  T  »  300°  and  60C^K  and 
10-percent  dissociation.  Here  one  sees  that,  at  both  30flP  and  6<K)Pk, 

Models  3  and  4  give  substantially  equivalent  results.  The  small  difference 
between  the  two  results  at  300°K  where  the  three-body  recombination  reaction 
dominates  is  due  to  the  more  exact  calculation  of  the  three-body  concentra¬ 
tion  that  is  possible  in  the  more  elaborate  Model  4.  There  is  no  signi¬ 
ficant  effect  of  the  atomic  nitrogen  regeneration  provided  by  one  of  the 
Type  (a)  secondary  reactions  included  in  Model  4. 


CAE  Report  No.  1054 
R-H22  Appendix  IV 

Volume  II 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Considering  the  simplest  case  of  recombination  in  pure  nitrogen, 
one  finds  that  the  process  is  controlled  by  a  three-body  reaction.  Its 
rate  thus  increases  with  pressure.  The  rate  coefficient  is  inversely 
proportional  to  temperature,  however,  so  that  the  lifetime  of  dissociated 
nitrogen  increases  with  increasing  temperature.  As  long  as  one  is  in¬ 
terested  in  the  asymptotic  time  behavior,  it  is  not  necessary  to  consider 
the  variation  in  the  third-body  concentration  with  time.  In  this  period 
n(N]'  time"*  and  is  independent  of  the  initial  atomic  nitrogen  concen¬ 
tration. 

2.  Under  typical  conditions  (STP,  10-percent  dissociation)  the 
asymptotic  time  period  is  reached  in  about  0.1  msec,  at  which  time  the 
atomic  nitrogen  has  decreased  by  two  orders  of  magnitude  from  Its  initial 
value. 


3.  The  effect  of  oxygen  on  the  nitrogen  recombination  depends  on 
temperature.  At  300°K  all  oxygen  reactions  are  unimportant  and  the  oxygen 
serves  simply  as  a  third-body  for  the  nitrogen  recombination.  In  this 
case  the  re.narks  In  (I)  and  (2)  above  apply.  By  600°K,  one  oxygen  reac¬ 
tion  has  become  sufficiently  important  to  completely  replace  the  three- 

*  body  recombination  as  the  most  important  atomic  nitrogen  loss  mechanism. 
The  reaction  is 

N  +  0  -  NO  +  0 

This  reaction  will  dominate  the  atomic  nitrogen  loss  in  the  asymptotic 
time  period  t  >  0. I  msec  and  it  will  result  in  a  much  greater  loss  rate 
than  would  be  the  case  for  three-body  recombination. 

4.  Although  detailed  calculation  of  the  energe.Ics  of  the  ammonia- 
ignition  process  has  not  been  undertaken,  it  would  ,iot  appear  feasible 

to  rely  completely  on  nitrogen  dissociation  for  this  function,  because  the 
major  part  of  the  recombination  will  occur  before  the  dissociated  nitrogen 
reaches  an  engine  cylinder. 

5.  This  work  does  suggest,  however,  that  RF  dissociation  of 
ammonia  may  be  of  Interest,  since  it  is  known  that  thermal  dissociation 
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of  ammonia  prior  to  ignition  is  •  ff-cti.o  i.i  improving  the  operation  of 
ammonia-fueled  engines.  The  RF  dlssoclator  ctin  be  expected  to  be  con¬ 
siderably  easier  to  control  than  a  thermal  dissociator,  and  should  not 
recuire  a  warm-up  period. 
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FIGURE  2.  HIGH- ENERGY  PULSER  DISCHARGING  IN  AIR 
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FIGURE  2.  HIGH- ENcRGY  PULSER  DISCHARGINGHN  AIR 

200(000  VOLT  S  '0  JOULE  S - 

R- !  122 


40 


CAE  Report  No.  10S-} 
Appendix  IV 
Volume  II 


FIGURE  3.  TWO  ELECTRODE  CONFIGURATIONS,  SHOWING  THE  TENDENCY 

TO  ARC  TO  THE  CYLINDER  HEAD.  THIS  WAS  SOLVED  BY  EXTENDING 
THE  ELECTRODES  DEEPER  INTO  THE  CHAMBER  AND  POLISHING 
THE  ROUGH  EDGES  OFF 
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FIGURE  8.  COMBUSTION  TIME  VERSUS  PRESSURE  PERCENT  OF  HYDROGEN 
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FIGURE  9.  COMBUSTION  TIME  VERSUS  YOTAL  PRESSURE  FOR  HYDROGEN 
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FIGURE  10.  COMBUSTION  TIME  VERSUS  DISCHARGE  ENERGY  FOR  HYDROGEN 
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FIGURE  II.  COMBUSTION  TIME  VERSUS  PRESSURE  PERCENT  AMMONIA 
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FIGURE  12.  COMBUSTION  TIME  VERSUS  TOTAL  PRESSURE  FOR  AMMONIA 


50 


8REAKFR  SYSTE 


CAE  Report  No.  1054 
Appendix  IV 
Volume  II 


TEMPE*IATJRE  [°ft] 


FIGURE  14.  COMBUSTION  TIMF  VERSUS  TEMPERATURE  FOR  AMMONIA 
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FIGURE  15.  PRESSURE  RATE  VERSUS  TEMPERATURE  FOR  AMMONIA 
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FIGURE  17.  MODEL  2,  10%  INITIAL  N9  DISSOCIATION 
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APPENDIX  IV- A 
Subroutine  PRECOR 

The  subroutine  PRECOR  used  for  some  of  the  numerical  calculations 
presented  here  is  particularly  well  suited  for  the  sets  of  coupled  non¬ 
linear  ordinary  differential  equations  encouncered  in  this  type  of  stuay. 
Because  of  its  general  utility,  the  following  discussion  of  the  subroutine 

and  a  listing  of  the  subroutine  is  presented  here.  This  material  is 

57 

abstracted  from  the  memorandum  of  Rank.  ' 
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INTRODUCTION 


k  c,3Urth-order  stable  predictor-corrector  technique  is  used  to  solve 
any  differential  equation  of  the  form 

dy  _ 


with  initial  conditions 


2  =  f  (x  ,  y) 
dx 


y(xQ)  -  yD 


THEORY 

In  order  to  simplify  what  follows,  y  and  f  will  be  taken  to  be 
scalars;  the  extension  to  N  dimensions  is  obvious.  Hamming  has  devised  a 
predictor-corrector  scheme  which  is  stable  and  has  the  advantage  that  only 
two  values  of  f  are  computed  at  each  step  (see  Ralston  and  Wilf,  *^and  Hamming**®) 

Using  the  notation  y  =  y  (x  )  ,  and  x  =  x  +  r,h  ,  the  equations 
•J  n  n  no  M 

for  Hamming's  method  are 


Predictor  P  =  y  ,  +  (2y  '  -  y  , '  +  2y  *) 

rs+l  n-3  3  n  n-l  'n-2 


Modi f ier 

r.+  l 

•B 

**n+’  * 

112  (p  . 
121  1  n 

Cn> 
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= 

f  (xn+i 
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Cn+I 
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An  estimate  of  the  truncation  error  is  given  r>y 


C  .  .  -  P  .  , 
ri+ 1  n+ 


-  1  —  I  ,y  fe  S 

-  36O  h  y  {§) 


x  .  <  §  <  x  , 

n-i  n+l 


Note  that  this  method  is  not  seif-starting,  for  it  requires  the  four 
previous  values  of  y  and  y' .  Therefore,  in  order  to  start  the  calcula¬ 
tion,  yx  .  y2  ,  and  yrj  must  be  computed  (y0  is  given).  A  fourth-order 
Runge-Kutta  technique  devised  by  A  Ralston^ is  used  to  compute  these 
va 1 ues . 

If  the  truncation  error  becomes  too  large,  the  interval  h  can  be 
halved  by  using  the  following  sixth-order  interpolation  formulas  for  the 
needed  values  of  y. 

V}  *  256  (8%  +  l35yn-l  +  ',0yn-2  +  yn-3> 

+  215  *  'Syn.2'> 


yn4  =  256  (l2yn  +  l35yn-l  +  l08yn-2  +  yn-3> 
+  (-3y  -  54y„  '  +  27y„.,') 


If  the  truncation  error  becomes  smaller  than  required,  the  interval 
may  be  doubled  if  enough  of  the  previous  values  of  y  and  y'  are  known. 


FCR1  Ft  AN  PROGRAM 


FLOW  CHART 

Figure  A-l  is  a  flow  chart  for  PRcCOR.  A  brief  description  of  each  sub¬ 
routine  which  appears  on  this  chart  follows: 

STnTuF  The  first  three  values  of  y  and  y'  needed  io  Start  the 

predictor-corrector  scheme  are  computed  by  a  Runge-Kutta 
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technique,  if  t>o  truncation  error  is  too  large,  the'  spaclnq 
h  is  halved  and  the  calculation  repeated. 


DIFIF 

This  subroutine  Is  called  by  STRTUP  to  compute 

of  the  truncation  error  after  the  first  three 

an  estimate 

calculations . 

CALC 

This  calculates  the  next  value  of  y  and  y', 
the  truncation  error. 

and  estimates 

HALF 

h  is  halved. 

DOUBLE 

h  is  doubled. 

ADJUST 

This  is  essentially  a  bookkeeping  subroutine; 

it:  rearranges 

the  values  of  y  and  y1,  to  prepare  for  the  next  calculation. 


The  control  variable  JJ  is  used  for  testing  to  find  out  if  there 
are  enough  previous  values  of  y  and  y'  stored  in  memory  in  order  to 
double  h. 

The  results  of  the  integration  ar  •  punched  on  cards  with  the  following 
format : 

X  1  y’  (x) 

2  y*M 


n 


The  estimate  of  the  truncation  error,  SUM,  is  oivei  by  a  weighted  sum 
of  the  absolute  value  of  the  diffcr<ncc  between  the  coinpo”r;nts  of  the  pre¬ 
dictor  and  corrector. 
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where 

n 

* 

the  dimension  of  y 

a. 

i 

•  = 

we i qhts 

pi 

= 

the  component  of 

the  predictor 

c' 

= 

the  i1*1  component  of 

the  corrector 

Two  conditions  are  provided  in  PRECOR  for  return  to  the  mainline 
program.  The  first  is  smallness  of  the  spacing,  relative  to  the  input 
parameter  HMIN;  singular  behavior  of  the  solution  will  usually,  be  the  cause 
of  this  exit.  The  other  exit  occurs  if  the  ratio  X/X..F  becomes  greater 
than  1,  where  XF  is  the  final  value  of  X  for  which  the  integration  is 
to  be  carried  out.  If  the  programmer  desires  other  criteria  for  exiting, 
e.g.  when  the  second  component  of  y  exceeds  5,  these  conditions  should 
be  Inserted  into  PRECOR,  and  the  existing  test  removed.  Note  that  in 
general  these  changes  will  affect  the  arguments  of  PRECOR,  but  will  not 
affect  the  COMMON  statement. 

Irrespective  of  the  truncation  error,  it  may  be  desirable  that  the 
spacing  h  never  exceed  a  certain  value.  This  value  is  given  by  HHAX. 


EXPLICIT  ARGUMENTS  OF  PRECOR 


X0 

Y0 

SMAX 

SMIN 

HMAX 

HMIN 

XF 

FUNC 


x  ,  initial  value  of  x 
o 

,y  ,  initial  value  of  y 

If  the  estimation  of  truncation  error  is  greater  than  SMAX, 
spacing<  is  halved. 

If  the  estimation  of  truncation  error  is  less  than  SMIN, 
spacing  is  doubled.  The  difference  between  SMAX  and  SMIN 
should  be  at  least  one  order  of  magnitude. 

Maximum  value  of  the  spacing  allowed 

Minimum  value  of  the  spacing  allowed 


:inal  value  of  x 

lummy  name  of  subroutine  to  compute  f(x,  y) ;  the  actual 
lame  used  in  the  argument  of  PRECOR  in  the  mainline  may  be 
iny.thing  except  FUNC,  PRECOR,  STRTUP,  BIFF,  CAL'C,  HALF, 
10UBLE,  dr  ADJUST.  This  name  must  appear  in  an  EXTERNAL 
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X 

W 

WO 

Cl 


statement  in  the  mainline  program  and  must  be  a  subroutine 
subprogram;  i.e.  a  CALL  statement  must  be  used  to  call  it. 
The  arguments  for  this  function  must  be  of  the  form  (X,  W, 
WD,  Cl),  where 

is  a  non-d imens ioned  variable; 

is  a  dimensioned  variable  with  the  same  dimension  as  Y  ; 

'ic-  a  dimensioned  variable  with  the  same  dimension  as  Y  I 

is  a  dimensioned  variable  whose  elements  may  be  used  as 

parameters  in  defining  f  ; 

and 


X,  W,  and  WO  must  be  related  by 


-WO  (I)  =  f  (X,  W{I)) 


IMPLICIT  ARGUMENTS  OF  PRECOR 
(CALLED  THROUGH  COMMON) 

H  h,  initial  value  of  the  spacing 

•  4 

N  Dimension  of  y  and  f  ,  i.e.  number  of  components 

A  Weights  needed  in  computing  SUM 

Cl  Parameters  required  to  define  f 

OTHER  REQUIREMENTS 

A  COMMON  statement  of  the  type  appearing  in  PRECOR  must  appear  in  the 
mainline  but  not  in  the  subroutine  used  for  calculating  r{x,  y)  . 

The  program  as  written  restricts  the  maximum  dimension  of  y  to 
eight,.  To  remove  this  restriction,  change  all  eight's  in  the  dimensioning 
information  to  the  dimension  of  y  required.  Do  not  change  the  seven's. 

This  program  must  be  used  with  the  Kingston  Fortran  li  compiler,  the 
restriction  is  imposed  because  the  I BM  1620  Fortran  |!  compiler  lacks  the 
EXTERNAL  and  BEAL  statements.,  A  l-lst-sng  of  the  program  follows.. 
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statement  in  the  mainline  program  and  must  be  a  subroutine 
subprogram;  i.e.  a  CALL  statement  must  be  used  to  call  it. 
The  arguments  for  this  function  must  be  of  the  form  (Xj  W, 
WO,  Cl),  where 

X 
V 
WD 
Cl 


X,  W,  and  WO 

.WD  (I)  -  V  (X,  W(I)) 


is  a  non-d imens ioned  variable; 

is  a  dimensioned  variable  with  the  same  dimension  as  Y  : 

*ic-  a  dimensioned  variable  with  the  same  dimension  as  Y  ; 

is  a  dimensioned  variable  whose  elements  may  be  used  as 

parameters  in  defining  f  ; 

and 

must  be  related  by 


IMPLICIT  ARGUMENTS  OF  PRECOR 
(CALLED  THROUGH  COMMON) 

H  h,  initial  value  of  the  spacing 

N  Dimension  of  y  and  f  ,  i.e.  number  of  components 

A  Weights  needed  in  computing  SUM 

Cl  Parameters  required  to  define  f 

OTHER  REQUIREMENTS 

A  COMMON  statement  of  the  type  appearing  in  PRECOR  must  appear  in  the 
mainline  but  not  in  the  subroutine  used  for  calculating  r{x,  y)  . 

The  program  as  written  restricts  the  maximum  dimension  of  y  to 
eight,.  To  remove  this  restriction,  change,  all  eight's  in  the  dimensioning 
information  to  the  dimension  of  y  required.  Do  not  change  the  seven's. 

This  program  must  be  used  with  the  Kingston  Fortran  li  compiler,  the 
restriction  is  imposed  because  the  f BM  1620  Fortran  jl  compiler  lacks  the 
EXTERNAL  and  REAL  s  tuteirr-rrts .  A  listing  of  the  program  follows,. 
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SUbROUT I Nb  PRhCOR ( XO . YO , SMAX ♦ SM I N . HMAX ,  HM 1 N  .  XF . FUNC 1 
COMMON  X»Y<7,8) »YD<7«8) .H.N.DtS) ,0S<8) ,A(8)  .5UM.JJ.C1  «50» 
DIMENSION  Y0<8) 

100  FORMAT  (  //7X  *<1  HX  •  1 2X  »  I  H  1  *  7X  .  4HY  t  I  )  /  ) 

99  FORMAT (£18.8* 16»fcl  7.8) 

98  FORMAT (  1 SX . I  6 .b 1 7.8 1 

PUNCH  1 00 
JJ=3 

CALL  STRTUP(XO.YO.FUNC.SMAX) 

X=XO 

DO  4  J= 1.4 

* 

1  =  1 

PUNCH  99. X »  I-.Y <  J+2, 1  ) 

IF  1 N— 1 14.4.1 

1  DO  2  I  *  2  «  N 

2  PUNCH  98.  I .Y( J+2. I  1 

4  X=X+H 

5  CALL  CALC  (FUND 
IF (SUM— SMAX 17.6.6 

6  CALL  HALF  .<  FUNC  1 

1 F (HM I N/H— 1 *0)8*5. 17 

7  1  =  1 

PUNCH  99.X. I »Y<  7. 1  h 
1F(N-1 1 10, 10,8 

8  DO  9  1=2^N 

9  PUNCH  98,  1  «,Y  <7,  1  1 

10  IF(XXXF-t.O)  11.17.17 

11  1F<JJ-1)  13.13.12 

12  JJ=JJ-1 
GO  TO  14 

13  IF { SUM— SM 1 N 1  15.15,14 

14  CALL  ADJUST 
GO  TID  5 

15  IF(r2.0*H/HMAX-l  .0)  16.16.14 

16  CALL  DOUBLE 
GO  TO  5 

17  RETURN 
END 
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STRTUP-  USES  RUNGE-KUTTA  METHOD  TO  COMPOTE  THE  F  I  RSI  A  VALUES  OF  Y 
SUBROUTINE  o  TP  T  OP  (  XO  .  YO  .FUNC  •  SMAX  )  ANP  F 

COMMON  X* Yt 7.B  > .YD( 7.B) .H.N.DtB  1  .DS(B 1 «A(8)»SUM« JJ.C1  <bU» 

DIMENSION  YO  (  8  )  *Wi'B)  «  WD  (  £)  .YD0<8)  .K 1  (8)  .  Kt'(  B  )  »K  J<  8  )  ,K.4(B> 

REAL  Kl  .K2.K.3.K4 

CALL  FUNC  (  XO  «  YO  •>  s‘DO . C  1  ) 

00  1  1  =  1  »N 

Y  (  3  «  I  ) = YO (  l ) 

1  YD <  3 «  I.)  =YDO (  1  ) 

7  X  =  XO 

DO  6  J=  1  *  J 
DO  2  1  =  1  »  N 
K.1  (  I  )  =  YDI  J4-2.  l  )  *K 

2  W( 1 )=Y< J4-2. 1 >+O.A*M t ti 
X=X4U.4*H 

CALL  FUNC ( X. W. WD. C 1 ) 

DO  3  I =1 vN 
K2( 1 )=H*WD< I ) 

3  W( 1 ) =W< I ) ~. 10J02239*K1 ( 1 )  4- .  1  b8  75964*K2< I > 

X=X+,0557372b*H 

CALL  FUNC (X.W.WD.Cl ) 

DO  A  1=1. N 
•K3< i )=H*WD( 1 > 

A  W<  1  )=W(  I  )  —»  U780772  1  *K  1  (  ],)-J.209724B*K2(  I  >4-3.8J2B648«K3C  1  > 
X=X4-.S442627b*H 
CALL  FUNCtX.W.WO.Cl > 

DO  5  1=1. N 

K4(T )=H*WD< 1 > 

5  W( 1 >=Y< J+2. 1 >4. 174 76028 *K1 ( 1 ) -.55 1 4Q066*K2 ( I )+ 1 .2055356«KJ C 1 )  +  . 171 
1 18478*KA<1 ) 

CALL  FUNC (X.W.WD.Cl ) 

DO  6  I =1 »N 
Y(U43. I >=W< I) 

6  YD( J+3. 1 ) =WD( 1 ) 

CALL  D1FF 

1 F ( SUM— SMAX )  8.9.9 
9  H=H«0.b 
CO  TO  7 

8  Rb  TURN 
bND 


DIFF—  COMPUTES  PRED.  COP.  D I FF ERENCES  AND  THEIR  WEIGHTED  BUM 
SUBROUTINE  DIFF 

COMMON  X.Y(7»8).YU(7»0).H.N.DlB).DS(8).A(8r. SUM . J J  « C 1  ( 50  > 

CONSl =0.962963 
C0NS2=0«375*H 
SUM=0*0 
DO  1  1=1. N 

DS  C  I  >=C0N5l*(  Y<6.  J  > ~Y .(,3...  I. )  — CONi>2* .(.Yt> ( 6..  .1  )  +J.O«y.D( b.  I  )43,.U*  YU  t  A  .  1  J  t 
lYO(3. I ) ) ) 

1  SUM=SUM4M(  1  >*At>SF  (DS(  I  )  ) 

RETURN 

END 
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C  CALC-  COMPUTES  NEW  VALUES  OF  Y  AND  F<X»Y> 

SUBROUTINE  CALC  <  FUND 

COMMON  X  «  Y  <  7  ♦  13  )  *YD<7.B>  .  H .  N«  D  (  8  >  ♦  DM  0  )  *  A  (  U  b»  SUM  *  J  J  4  C  1  C£jO) 

D  1  ML  NS  I  ON  M<8)  *MD<8)  «P<8)  *C<B>»WCB>  *tfD<0.> 

REAL  M*MD 

CONS l = i • J3333J J*H 

coNS2=*92b6J98J 

C0NS3= 1 *0— C0NS2 

SUM=0«0  H 

DO  1  1=1 iN 

PC  1  )=Y<3* I 1+CONSl* <2oO*YD<6* 1 >-YDCb« 1 >+2.0*YD<4 » t ) ) 

1  MM  )=P<  1  )-C0NS2*DS<  1  } 

CALL  f*  UNC  <  X*M«MD*C1  } 

DO  2  I,=  l  «N 

CC  1  )=0«125*(9.0ftY<6«  1  )-YU«  1  )+3.0»H#(MD<  1  >+2*U«YD(6«  1  >-YD<b»  1  )  >  > 
DM  )  =P(  1  I  -C  <  I  ) 

SUM=SUM+A<  I  )»aBSFO( 

2  W< i )=C( 1 )+C0NS3*U< t) 

CALL  FUNC ( X* W*  WO  *C 1 ) 

DO  3  I  =  I  .  N 

Y  (  7  « I  ) =W<  1  ) 

3  YD<  7. 1  )=WUl  l  ) 

RETURN 

END 


C  HALF-  HALVtS  THE  SPACING 
SUBROUTINE  HALF! “UNO, 

COMMON  X  *-Y  (7i8)  *  YDC  7*8>.»H«N*D<8  )  «DS(8)  iA(BI  »SUM.JJ.C1  CbO) 

DIMENSION  W(8).WU<8) 

CONS 1  =  • 0039062b 
DO  1  1=1.N 

YA=CONSl* C80.U*Y<6. 1  )  +  i3b*0*YC5»  l  >+40.0»Y<4*  I  )+YC3*  1  <>+H*<-lbiD«YuC 
16*1  ) +90  »0*YD ( b  *  I  )+lb*0#YD<4» 1  J)  > 

YB=CONSl*<  12*  0*Y ( 6  *  I)  +  1 Jb • O* Y ( b . I  )  +  1 08»0* Y { 4 ♦ I  )  +  Y ( J* I  ) +H* ( — J *u *  YU < 
16* I )-b4*U*YD(b* 1 )+27.0#Y0<4* 111) 

YD (2* I >=YO(4, 1 ) 

YD ( 4  » I  >=YDib« I  1 
Y(2s  l  );=  Y  (  4  »  l  ) 

YC4  «ij  )=Y(b«  I  ) 

YiS. 1 >=YA 
1  YC3.1J=YB 
H=0.5*H 
FA=3.0 

DO  3  J=bi3«~2 
DO  2  1=5 «N 
*  W(  i  1.-.Y5  J*  I  ) 

XC=X-FA*H 

CALL  FUNC  (  XC «  W  ■>  WD  *  C 1  ) 

T'A‘-yiw 
DO  3  1=1)N 


{cont ’d] 
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3  YD< J. 1 > =WD< I » 
X  =  X-H 
JJ =2 
re  turn 

END 


C  DOUBLE-  DOUBLES  THE  SPACING 
SUBROUTINE  DOUBLE 

COMMON  X* Y<  7*8 ) » YD<  7*8  >  *H«NiU<8  >  «DS(U  > ♦ A18)  ♦  SUM  *  JJ  *  Cl  t  50  ) 
DO  1  1=1 «N 

Y<6* I )=Y( 7. 1 ) 

Y< At  I )=Y<3* I ) 

Y  <  3  *  IJ=Y<  1  *  1  ) 

YD<  6* I )  =YD<7« 1  1 
YD ( 4  4  I  }=YD<3« I  J 
YD  1 3  *  1  >=YD( l ♦ 1  > 

1  DS< 1 )=D( 1 ) 

H=2.0*H 
J  J  =  3 
X=X+H 
RETURN 
'  END 


C  ADJUST-  rearranges  Y  AND  F<X«Y>  ARRAYS  TO  PREPARE  FOR  NEXT  PASS 
SUBROUTINE  ADJUST 

COMMON  X«Y17«8) *YD«7*8) .  H.N.D  <  B  )  *DS  (  8  >  »  A  ( 8  >, «  SUM  *  J  J  ,  C  H  50  > 

DO  1  1  =  1  *  N 

DS< 1 )*D< l ) 

DO  1  J=JJ*6 
YOU.  1  >=YD< J+l  «  1  1 
1  Y< J» l »=Y( J+l . I ) 

X=X+H 

RETURN 

END 
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ABSTRACT 

This  study  Is  to  determine  the  feasibility  of  obtaining 
hydrogen  for  Improving  combustion  of  ammonia  In  Internal  combustion 
engines. 

important  recorrWngttfon  reeet’ons  and  their  rate  constants 
were  obtained  from  the  literature,  end  the  combined  rate  equations 
solved  numerically  with  ©  digital  computer* 

Stefidy~3i«te  concentrations  of  hydrogen  and  hydrazine  were 
found  with  Initial  dlssssoctetion  from  1%  to  10%  and  temperatures 
from  0°  to  I00°C. 

Yields  were  found  to  be  too  low  for  practical  use  In  an 
ammonia  engine. 


Keywords:  Engine,  internal  combustion 

Ammonia 

Fuels,  internal  combustion 
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STEVENS  INSTITUTE  OF  TECHNOLOGY 

DAVIDSON  LABORATORY 
CASTLE  POINT  STATION 
HODOKEN,  NE\W  JERSEY 


AMMONIA  DISSOCIATION  BY  RADIO-FREQUENCY  ENERGY 

INTRODUCTION 

I  **' 


Background 

Within  the  framework  of  a  program  to  determine  the  suitability  of 
ammonia  as  an  engine  fuel,  CAE  conducted,  along  with  other  work,  spark- 
ignition  engine  tests  in  which  a  thermal  dissociator  produces  free 
hydrogen.  Stevens  Institute  has  investigated  the  use  oF  high-energy 
spark  discharge,  and  radio-frequency  ionization  of  the  nitrogen  in 
the  intake  air,  as  means  of  increasing  ignition  energy. 

Experience  with  the  thermal  dissociator  has  shown  that  operation  is 
satisfactory  when  sufficient  hydrogen  is  generated.  But  during  starting 
and  certain  transient  conditions  the  hydrogen  produced  is  insufficient 
for  operation.  Also,  there  is  evidence  of  limitations  on  catalyst 
1 ifetime. 

Therefore,  a  look  at  means  other  than  thermal,  for  the  dissociation 
of  arrmonia,  -is  in  order.  The  primary  concern  of  this  work  is  to 
determine  whether  an  RF  dissociator  can  substitute  for  or  operate  in 
conjunction  with  the  current  thermal  unit,  to  alleviate  present  shortcomings. 

Character  of  Work 

Although  in  the  final  analysis  the  merits  of  this  scheme  must  rest 
upon  an  experimental  demonstration,  analytical  studies  at  the  first  stage 
will  be  useful  in  establishing  the  feasibility  of  the  approach,  and 
helpful  in  defining  the  general  nature  of  arty  experimental  work  to  be 
undertaken. 

The  purpose  of  this  work  is  to  establish,  analytically.,  the 
bility  cf  ammonia  dissociation  by  the  use  of  RF  or  other  energy  and  to 
estimate  the  size  and  power  requirements  of  a  suitable  device. 


1 


CAE  Report  No.  1054 
Appendix  V 
Volume  II 


R-1210 


ANALYSIS 


In  recent  years,  radio-frequency  discharges  have  beer,  used  as  ion 
sources  in  plasma  physics,  as  high-temperature  torches  for  the  welding  of 
refractory  materials,  and  as  detectors  of  ionizing  radiations,' 

The  oscillatory  field  may  be  produced  in  the  gas  by  electrodes  to 
which  the  high-frequency  potential  is  applied,  or  by  induction  from  e 
coil  carrying  an  oscillatory  current*  or  the  gas  could  be  contained  in  a 
wave  guide  or  resonant  cavity. 

The  mechanism  of  radio-frequency  dissociation  operates,  in  general, 
in  the  following  way.  Initially,  a  small  amount  of  the  oas  is  ionized, 
producing  free  electrons;  these  electrons  are  accelerated  by  the  oscilla¬ 
ting  electromagnetic  field  and  collide  with  the  molecules  of  the  gas;  the 
energy  of  collision,  ,f  sufficient,  will  cause  dissociation  of  the  molecules 
of  the  gas;  and  the  dissociated  products  will  then  seek  to  recombine,  the 
eventual  state  of  the  gas  being  dependent  on  the  kinds  of  recombination 
reactions  possible,  on  the  relative  probabilities  of  the  recombinations, 
and  on  the  stability  of  the  recombination  products. 

The  amount  of  dissociation  i3  dependent  on  the  RF  power  available, 
the  efficiency  with  which  the  RF  is  couplod  into  the  gas,  and  the  above 
recombination  considerations-  The  coupling  efficiency  depends  on  among 
other  things,  the  frequw'y  «f  the  RF  source.  If  the  electrons  are 
accelerated  to  their  ma,  .muni  velocity  ir»  about  one  mean  free  path,  their 
coll£<ional  energy  will  be  highest.  Unfortunately,  at  atmospheric  pressure, 
this  corresponds  to  an  impossibly  high  frequency. 

The  primary  objective  of  this  study  is  to  establish  the  yield  of 
free  hydrogen  in  the  mixture,  as  a  function  of  the  temperature  and  the 
density  of  the  reactants.  From  the  system  of  recombination  reactions  that 
can  occur,  sets  of  coupled  ordinary  non-linaar  differential  equations  are 
constructed?  these  govern  the  time  and  tesperaturu  dependence  of  the  con¬ 
centrations  of  ths  various  Churiltal  speoiss.  Mie  equations  win  then  be 
solved  by  the  use  of  a  high-speed  digital  computer,  to  determine  what 
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environmental  and  initial  conditions  must  be  satisfied  if  an  adequate 
concentration  of  hydrogen  is  to  be  maintained  in  the  emmonia-air  mixture. 

In  addition  to  hydrogen,  hydrazine  (^H^)  evo^ves*  Although  not 
of  primary  interest,  the  hydrazine  concentration  will  also  be  examined, 
to  determine  whether  or  not  the  effects  of  this  component  are  significant. 
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PROCEDURE 

The  work  of  this  project  was  divided  into  the  five  major  phases 
1 isted  below: 

(1)  Assemble  data  on  the  reaction  kinetics  and  rate  coefficients 
of  the  ammonia  system.  The  literature  search  on  these 
kinetics  was  performed  by  Or.  frank  Jones. 

(2)  Set  up  mathematical  models  of  the  reaction  system,  based  on 
the  rate  equations  obtained  in  phase  (1)  above. 

(3)  Solve  the  sets  of  equations  on  a  high-speed  digital  computer, 

2 

using  subroutine  PRECOR. 

(4)  Interpret  the  results  as  they  apply  to  the  physical  problem. 

(5)  Perform  a  check  on  the  accuracy  with  which  the  mathematical 
analysis  describes  the  physical  system 
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RESULTS 


Reaction  Mechanism 

Tht  collision  between  ammonia  molecules  and  electrons  accelerated 

in  an  electric  field  produces  excited  ammonia  molecules  which  may 

ultimately  decompose  when  the  excitation  energy  is  sufficiently  high. 

3 

Anderson  et  al  have  shown  that  the  energy  required  for  dissociation  in 
an  RF  discharge  is  very  ciose  to  tie  N-H  bond  energy  ir.  NH,  (4.5  eV), 
implying  that  the  primary  products  and  H  and  NIL  radicals  in  their  ground 
states}  thus 

nh3 .  m2  +  H  (0) 

Other  decomposition  modes  such  as  IlH^  —  I1H  +  require  higher  energies 

and  can  therefore  occur  only  to  an  insignificant  extent.  Reaction  (0) 

is  also  the  sole  .node  of  decompos i tion  in  the  photochemistry  of  ammonia 

0  4.5 

for  light  of  wavelengths  above  1500  A.  ’  . 

The  primary  decomposition  products.  H  and  tili^.  may  interreact  in 
three  different  ways.’ 

H  +  H  +  M  — —  H2  +  11  (1) 

H  +  +  M  — -  +  H  (2) 

2  3 

NH2  +  MH2  +  M -  N2H4  +  H  (3) 


where  H  is  any  third  body. 


Reaction  (1)  is  exothermic  by  an  amount  equal  to  the  H-H  bond  energy,  so. 
the  H2  molecules  thus  formed  have  a  high  probability  of  splitting  up 
again  unless  the  energy  can  be  dissipated  to  other  molecules  in  the  system. 


The  rate  of  reaction  (1)  therefore  depends  not  only  on  the  concentration 


nf  M 


21 


Iso  on  the  tvtSi  gos  pressure* 


o 

n  —  t. 

■uut;  ”  \r*/« 


Similar  arguments  hold  for  reactions  (2)  and  (3h  however,  in  these  two 
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cases  the  respective  products  NH^  and  K2H^  are  complex  enough  to 
allow  for  storage  of  some  of  the  released  energy,  in  other  bonds  within 
the  molecule  -  hence  the  molecule  can  exist  fer  a  finite  length  of  time 
before  decomposing  .  The  rates  of  these  reactions  should  depend  on 
total  pressure  only  at  low  pressures  and  should  reach  a  limiting  value 
as  the  pressure  is  increased.  Both  reactions  (2)  and  (3)  are  therefore 
independent  of  total  gas  pressure  at  pressures  above  200-400  rrm  Hg.  For 
1  atm.  the  rates  are  given  in  Table  !a»^ 


If  a  third  body  is  not  present  when  two  UH-  radicals  collide*  the 

7  g  i  g 

reaction  below  mty  occur,'  with  a  rate  constant  k,  -  0.46  x  10  Hole 
a  3a 

(cf.  k.3  =  2.5  x  10.). 


sec 


NH2  +  NH2  -----  NHj  +  NH  (3a) 


g 

However,  the  reaction  of  HH  with  NH^  to  form  hydrazine  is  fast,  so  that 
the  net  result  is  stoichiometrical ly  equivalent  to  reaction  (3). 


NH  +  NH3  — —  Nj,H^ 


The  temperature  dependence  of  reaction-rate  constants  generally  obeys 
the  Arrhenius  equation  k  =  A  where  k  is  the  rate  constant  at 

obsolute  temperature  T,  E  is  the  Arrhenius  activation  energy,  and  A  is 
a  temperature-independent  pre-exponential  factor.  Reactions  between  free 
radicals  are  always  found  to  have  very  low  or  zero  activation  energies, 
and  therefore  the  rate  constants  for  reactions  (1),  (2) ,  and  (3)  arc  In¬ 
dependent  of  temperature. 

H  and  MH^  radicals  undergo  many  coMisic.,s  with  ammonia  molecules 
at  1  atm.  The  reaction  of  NH2  with  anmonia, 

nh2  +  nh3 -  kh3  *  nh2 

is  thermoneutral  and  can  be  disregarded  since  it  does  not  lead  to  chemical 
change.  The  reaction  of  H  with  ammonia  ha?  hn»n  the  subject  of  s  greet 
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deal  of  discussion?  it  is 


H  +  NH3  — —  H2  +  NH2  (4) 

Oixon^  and  Taylor  and  Jungers^  found  no  evidence  for  reaction  (4)  over 
the  temperature  range  25  -  260°C  and  placed  a  lower  limit  on  the 

activation  energy  of  E  s  8.5  k-cal/mole.  But  Jones  in  the  radiation 

12  d  13 

chamistry  and  Takamuku  and  Back  in  the  mercury-sensitized  photolysis 

explained  increased  yields  at  temperatures  up  to  350°C  on  the  basis  of 

reaction  {4).  Using  H  atoms  generated  in  a  hi gh-frecyjenry  discharge  in 

H-,  Volpi*4  estimated  an  activation  energy  of  10  -  15  k-cal/mole  and 

*  1213  15 

pre-exponential  factor  of  between  10  and  10  5  cc/mole  secj  Oganesyan  3 

determined  an  activation  energy  of  13*7  k-cal/mole  but  did  not  measure 

the  pre-exponential  factor.  Probably  the  best  determinations  were  made 

1 6 

by  Farkas  and  Melville  for  the  02-NH3  exchange  reaction  (Ea“  11-1 
k-cal/mole  and  A  *  2  x  10^  cc/mole-sec).  In  contrast  to  these  values, 
very  low  activation  energies'  of  1-2  k-cal/mole  have  often  been  erroneously 
assigned  to  reaction  (4),  but  these  values  must  refer  to  the  over-all 
process  of  hydrogen-formation  In  anmonia  and  not  to  reaction  (4). 

H2  *  NH2 - H  +  NH3  (4a) 


If  reaction. (4)  is  important,  then  (4a),  the  reverse  reaction,  must 
also  be  considered.  The  equilibrium  constant  ■‘or  the  exchange  of  the 
hydrogen  atom  between  and  as  measured  for  example  with  the 


deuterium  isotope  0,,  + 


HO  +  WLD,  has  been  determined  to  be 
19  L 


J  19  i 

K  =  3»0,  independent  of  temperature  .  Since  K  =  the  rate  constant 

for  (4a)  is  K,  1 


'4a  =  \  k4* 


Secondary  Reactions 

Hydrazine  formed  in  (3)  and  (3a)  may  be  destroyed  by  H  and  NH^ 
radicals? 
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H  +  N2H4 -  H2  +  N2H3  (5) 

nh2  +  h2h4  — —  nh3  +  n2h3  (6) 

14 

The  rate  constant  for  (5)  Has  been  measured  by  Volpi,  but  it  has  not 
yet  been  possible  to  measure  the  rate  constant  for  (6),  although  it  is 
expected  that  kgfltkg,  However,  the  results  indicate  that  (6)  can  be 
neglected  in  pure  arnnonia,  because  of  the  rapid  disappearance  of  NH2 
radicals  by  reaction  (3). 

The  unstable  N_H_  radica.  formed  in  (5)  or  (6)  may  be  further 

^  $  j  2 

degraded  by  reaction  sequence  (A), 

H  +  N2H3  — —  H2  +  N2H2 

h  +  n2h2 -  H2  +  N2H 

h+n2h  - —  h2  +  n2 

14  9 

although  Sequence  B  has  also  been  suggested,  "with  k„  >  3»0  x  10 

•  1  m\  “ 

Mole*1  sec  ,  so  that 

»2H3  +  N2H3  —  N2  +  2NH3 

Sequence  B  must  itself  be  the  net  result  of  a  series  of  steps  for  which 
k0  is  the  rate  constant  for  the  slowest  reaction.  In  either  case, 
reaction  (5)  is  the  rate-determining  step  whether  N2H3  decomposes  by 
A  or  B,  and  much  greater  than  or  Kg j  but  the  rate  of  H2  production 

by  the  sequence  (5)~(A)  will  be  four  times  that  for  the  sequence  (5)-(B). 

Kinetic  Rate  Equations 

The  reaction  mechanism  discussed  above  is  summarized  in  Table  1, 
and  the  appropriate  rate  constants  are  compiled  in  Table  ta.  The 
appropriate  differential -rate  equations  were  formulated  and  are  presented 

in  Tahlfl  ?. 


0 
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These  coupled  differential  equations,  with  initial  conditions 
specifying  the  level  of  first  dissociation  products,  were  then  solved 
on  a  Univac  1105  computer  with  subroutine  PRECOR. 

An  assumption  inherent  ir,  the  method  is  that  the  dissociation 
reaction  is  separated  from  the  recombination  reactions-that  is,  no  effort 
has  been  made  ito  examine  the  effects  of  simultaneous  dissociation  and 
t ecombi nation.  It  is  felt  that  this  assumption  is  of  negligible  consequence. 

At  first  glance,  the  reaction  mechanisms  would  lead  one  to  believe 
chat  the  Model  A  system  would  yield  a  considerably  greater  quantity  of 
hydrogen  then  Model  B.  Actual  computation,  however,  shows  that  the 
difference  -'n  hydrogen  yield  between  the  systems  is  insignificant". 

Since  the  reactions  are  temperature  dependent,  solutions  were  ob¬ 
tained  at  0°,  20°,  50°,  and  100°  C  (32°,  68°,  122°  and  212°F). 

It  was  not  feasible  to  solve  for  initial  dissociation  concentration 
in  terms  of  the  final  products,  and  it  was  therefore  necessary  to  assume 
a  set  of  initial  conditions  such  that  the  range  of  the  resulting  steady- 
state  hydrogen  concentration  spanned  the  region  of  interest. 

Steady-state  concentration  of  Hg,  and  is  presented  in 

terms  of  gram-moles  per  liter  in  Table  3»  mole  (volume)  percent  in  Table 
4,  and  weight  percent  in  Table  5*  The  column  labeled  is  explained  in 
the  section  entitled  "Oiscussion".  ~ 

Figure  1  shows  typical  concentration  vs.  time  behavior  on  a  linear 
scale.  This  is  shewn  for  reference  only,  since  the  crowding  on  the  left 
of  the  diagram  does  not  permit  accurate  observation  of  differences  in 
behavior. 

A  morn  useful  presentation  is  a  log-log  plot  of  concentration 
against  titfse.,  ac  in  Figures  2  through  13* 
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PRECISION 


In  order  to  correlate  the  analytical  study  with  the  a«~  ual  physical 

phenomena,  and  to  ensure  that  the  mechanism  selected  is  in  fact,  the  one 

followei  by  the  reaction  in  this  temperature-density  domain,  the  initial 

IP 

conditions  and  temperatures  employed  in  an  experimental  study  '  of 
ammonia  dissociation  by  radiation  were  inserted  in  the  computer  program. 

At  temperatures  of  23°  and  300°C,  with  initial  concentrations  of 

-3  -6 

H  and  NH^  of  3  x  10  percent ,  Jones  obtained  1,0  and  1.9  x  10  gram 

moles  per  liter  of  respectively,  and  an  increase  in  hydrogen  con¬ 

centration  by  a  factor  of  3*7* 

m 

Analytically,  at  those  same  points,  we  obtained  0.65  and  0.86  x  10 
gram  moles  per  liter  of  and  an  increase  in  H2  by  3*79- 
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DISCUSSION 


References  to  Tables  3  through  5  show  that  hydrazine  (NgH^) 
production  is  quite  significant.  In  fact,  the  reaction  efficiency  with 
respect  to  hydrogen  lies  in  the  vicinity  of  10%,  and  about  16%  for 
hydrazine. 

Since  no  conclusive  data  on  the  effect  in  engine  operation  of  a 
mixture  of  ammonia,  hydrogen,  and  hydrazine  is  available,  it  is  not 
possible  to  make  very  strong  direct  statements  regarding  such  a  mixture. 

The  steady-state  concentrations  calculated  in  Tables  3  through  5 
are  based  on  a  homogeneous  gas-phase  mechanism.  Heterogeneous  and  surface- 
catalyzed  reactions  have  not  been  consdiered.  Although  H 2  is  a  stable 
product  and  is  not  expected  to  be  influenced  by  metal  surfaces.  l^Hj,  is 
known  to  be  readily  decomposed  by  metals  even  at  low  temperatures.  The 
thermo-  dynamically  favored  mode  for  hydrazine  decomposition  on  a  surface 
is 

3H2H4 - "  N2  +  4NH3 

but  this  stoichiometry  is  achieved  only  in  aqueous  solutions  or  on  quartz 
surfaces  at  high  temperatures.  On  metal  surfaces,  is  usually  pro¬ 
duced  in  addition  to  Ng  and  There  is  no  agreement  on  the  stoichi¬ 

ometry  of  the  metal -catalyzed  decomposition  in  the  gas  phase,  and  results 

on  the  same  surface  are  often  found  to  vary  with  the  number  of  times  the 

20 

experiment  has  been  carried  out.  The  maximum  yield  from  the  de¬ 
composition  of  hydrazine  corresponds  to 

N2H4 .  N2  +  2H2 

and  the  yields  in  column  5  of  Table  5  were  calculated  on  this  basis. 

A  straightforward  calculation  proceeding  from  the  eV  first 
dissociation  energy  of  ammonia  yields  the  information  that  3*21  kilowatts 
are  necessary  to  dissociate  on  pound  per  hour  of  ammonia. 
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Assuming  a  brake  specific  arnonia  consumption  of  0.75  1b/8hp  hr.  a 

40-hp  engine  would  consume  30  1b/hr.  If  10%  dissociation  were  sufficient. 

(h!  =  0.56%  by  weight)  then  we  wouid  require 
2 

it 

33  NH3 .  30  NH3  +  3  NH3 

3  tf/hr.  .3.21  Kw  W/H  ~  9.65  Kw 

coupled  to  the  gas. 

*  22 

Cobine  and  Wilbur  of  G.E.  Research  constructed  an  RF  torch  and 
found  efficiencies  in  coupling  to  polyatomic  gases  of  70%.  nearly  all 
the  energy  going  into  dissociation. 

Thus,  an  over-all  efficiency  cf  55%  should  not  be  difficult  to 
realize. 

The  power  requirement  of  17-6  Kw  leads  to  the  untenable  result  that 
approximately  50%  of  tfte  total  engine  output  is  required  to  generate  the 
hydrogen. 

As  for  physical  size,  the  1  Kw  torch  of  Cobin  and  Wilbut  had  an 
rf  generator  8  5/8  inches  in  diameter  and  a  wave-guide-coupler  20-inches 
long. 


12 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Due  to  the  very  small  yield,  this  scheme  would  not  be  practical  as 
the  sole  source  of  hydropen  for  an  engine.  It  is  possible  that  this  scheme 
may  find  use  as  a  starting  and  warm-up  aid,  but  the  complexity  of  a  source 
of  high  frequency  is  discouraging. 

Therefore,  it  Is  recommended  that  this  approach  be  abendonned. 
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TABLE  1 

System  of  Chemical  Reaction,  Which  Govern  the  Recombination  Ij 


0, 

NH3  NH2  +  H 

H 

}. 

H  +  H  +  M  --->•  H2  +  M 

11 

2. 

H  +  NH0  +  H  — >  NH,  +  H 

U 

3. 

NH2  +  NH2  +  M  N2H4  +  H 

U 

3a. 

NH?  +  N!H2  —  MH3  +  NH  ; 

NH  +  NH,  — *»  N  H, 

3  2  4 

0 

4. 

H  +  NH3  — H2  +  NH? 

r,| 

u 

4a. 

H2  +  NH2  — H  +  NH3 

5. 

H  +  — *-  H2  +  N2H3 

!l 

A 

B 

j  f 
1  > 

H  +  ^2^3  — 

"  ^  H2  +  N2H2 

N2H3  +  N2»3  — 

2  NH3  +  N2 

,  i 

ij 

H  +  N2H2  — 

— — ■?»  H2  +  NjH 

ii 

H  +  N„H  — 

C 

-  -  -  H?  +  N2 

a 

3 
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TABLE  la 

Rate  Constants  of  the  Reaction?,  In  Table  1 


1.  1 3 .  +2.  x  103  mole2  sec”’  (room  temperature)  (Ea 

2.  12.  -1.  x  103  mole*  sec"*  (room  temperature  limit) 

3.  2.5  -.2  x  103  mole*  sec  *  (room  temperature  limit) 

a  -1  _i 

3a.  0.46  x  10"''  mote  sec 

4.  ,„I0  „-l0.7k-d/RT  mole*  sec"1 


4a.  1  . 

3  k4 


_  .  _  1n8  -2.0  k-cal/RT  ,-l 

5.  j.5  x  10  e  mole  sec 


-1 


A.  3-5 


,.8  -2.0  k-cal/RT 

10  e 

9  -1  -1 

!0  mole  sec 


-1  -1 
mole  sec 


0) 

(Ea  =  0) 
(Ea  *  0) 
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Condi tion 


0° 

1% 

4.415 

X 

io“2 

3.8 

X 

ic’5 

7.2 

X 

10  * 

5% 

4.237 

X 

io'2 

2.25 

X 

to-11 

3.75 

X 

~l± 
10  4 

10% 

4.014 

X 

JO*2 

3.85 

X 

io'*1 

6.99 

X 

O 

4 

20° 

l% 

4.118 

X 

ro*2 

4.0 

X 

IO-5 

6.8 

X 

O 

a 

vn 

5% 

3-952 

X 

io"2 

J:.0 

X 

,o'4 

3-4 

X 

IO-4 

10% 

3-744 

X 

10"2 

4.J 

X 

io-" 

6-6 

♦ 

IO-4 

O 

o 

i% 

3-732 

X 

io“2 

3.80 

X 

io'5 

6-7 

X 

10  3 

5% 

3-581 

X 

10’2 

2.0 

X 

io-4 

3-1 

X 

io"4 

10% 

3-393 

X 

10"2 

3-95 

X 

-4 

10  4 

6.0 

X 

IO**4 

o 

o 

o 

1% 

3-237 

X 

io-2 

4.0 

X 

-5 

10  5 

7-5 

X 

10'5 

5% 

3-106 

X 

10-2 

1-7 

X 

io'4 

2.8 

X 

io'4 

10% 

2-943 

X 

io”2 

3-78 

X 

10-4 

5-4 

X 

io*4 
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TABLE  4 

Concentration  In  Volume 

Percent 

Condition 

h 

V* 

0° 

1% 

0.085 

0.161 

5% 

0.505 

0.84 

10% 

0.865 

1.56 

20° 

1% 

0.096 

0.163 

5% 

0.48 

0-82 

10% 

0/98 

1-58 

50° 

1% 

0.101 

0.18 

5% 

0-53 

0-82 

10% 

1.05 

1-59 

100° 

1% 

0.122 

0.23 

5% 

0.52 

0.86 

10% 

MS 

1.65 

-20- 
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TABLE  5 


Weight  in  Percent 

Relative  to  Ammonia 

Condi tion 

X  -  h2 

x  - 

i 

o 

o 

0.0101 

0-306 

5% 

0.0625 

1.66 

10% 

0.  M3 

5-27 

0.52 

20°  1% 

0.0115 

0.31 

5% 

0.0535 

1.62 

10% 

'•1.29 

3. 31 

0.54 

50°  1% 

6  ,2 

0.337 

5% 

0.065 

1.625 

10% 

0.137 

3-32 

0.554 

100°  1% 

0-0145 

0.436 

5% 

0.0645 

1.695 

10% 

0-152 

3-^5 

0.584 

wt  x 
"  wt  MHj 


NOTE:  % 
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